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Abstract of J P 1 1 1 08845 (A) 

PROBLEM TO BE SOLVED: To obtain information 
on a resolution and a composition in a depthwise 
direction of an object to be inspected, by projecting 
a waveform including electromagnetic radiation 
pulses in a 10<12> Hz frequency range and 
measuring a related delay time of pulses reflected at 
the object SOLUTION: Electromagnetic radiation T- 
ray pulse beams in a T-ray frequency range 
generated by a tera Hz (T-ray) wave transmitter 103 
are reflected at a first reflecting surface 1 1 1 and a 
second reflecting surface 1 1 3, brought in an 
inspection object 107 and reflected. The reflected 
beams are reflected at a third reflecting surface 109 
and received by a T-ray wave detector 105. In this 
case, a reflection waveform includes a plurality of 
distorted pulses with time of various amplitudes and 
polarities because of the reflection at a 
discontinuous dielectric body in the inspection object 
107. An arrival time of each of the plurality of 
detected reflected pulses is measured at an 
operation process. A position of an interface of the 
dielectric body in the inspection object 107 in a 
propagation direction is determined from related 
delay times. Information on a resolution, a 
composition in a depthwise direction of the 
inspection object 107 is obtained in this manner. 
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(54) Method and apparatus for terahertz imaging 

(57) Certain material and objects can be character- 
ized by a frequency-dependent absorption, dispersion, 
and reflection of terahertz transients in signals which 
pass illuminate the material or object. The present ter- 
ahertz imaging system analyses that frequency depend- 
ence in the time-domain by collecting that transmitted 



signal propagating through the object and then process- 
ing the information contained in those signals for every 
point or "pixel" on that object. This is a non-invasive im- 
aging technique that is capable of differentiating be- 
tween different materials, chemical compositions, or en- 
vironments. 
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Description 

This is a division of application Serial No. 
08/388,933 filed February 15, 1995 entitled "Method 
and Apparatus for Terahertz Imaging". 5 

Technical Field 

This invention relates to spectroscopy in the tera- 
hertz frequency range and, more particularly, to a meth- 10 
od and apparatus for creating images of objects with sig- 
nals in this frequency range. 

Background of the Invention 

Terahertz time-domain spectroscopy ("THz-TDS") 
is a very powerful spectroscopic technique in the far- 
infrared spectral region. Terahertz radiation has been 
generated and detected using optically gated transmit- 
ters and receivers such as photoconducting dipole an- 
tennae as described in P. Smith et aL, IEEE J. of Quan- 
tum Electronics, Vol. 24, No. 2, pp. 255-260 (1 988) and 
N. Katzenellenbogen et aL, Appi. Phvs. Lett., Vol. 58, 
No. 3, pp. 222-224 (1991). With these techniques, ter- 
ahertz spectroscopy offers a reasonably good signal-to- 
noise ratio ( up to approximately 10 4 ); it can be per- 
formed without special thermal stabilization apparatus 
such as cooled detectors; it can be realized in a compact 
system; and it offers a transmitter and detector technol- 
ogy which is compatible with integrated circuit technol- 
ogy. 

Numerous experiments using terahertz time do- 
main spectroscopy have been performed on solids, liq- 
uids, and gases. Some experiments have analyzed the 
spectrum of a terahertz signal affected by carriers in 
semiconductors and superconductors. Other experi- 
ments have performed terahertz time domain spectros- 
copy on water vapor as well as N a O gas. Still other ex- 
periments have reported terahertz time domain spec- 
troscopy of chemical compounds in the liquid phase. In 
all these experiments, the terahertz signal was transmit- 
ted through the object under study in a single illuminated 
volume region (usually 25 mm in diameter) to provide 
the spectral information about that homogeneous re- 
gion. 

Summary of the Invention 

I have recognized that the time domain spectrosco- 
py and, more particularly, terahertz signals can be used 
for imaging objects by collecting individual signals prop- 
agating through distinct (spatially separate) points on 
the object and processing these signals to create the 
image of the object. It is also possible to focus the signal 
source on the object at distinct points and scan the 
source and detector in synchronism across the object 
in a pattern transverse to the propagation direction. Ad- 
ditionally, it is possible to cause to source to bathe the 



entire object with substantially parallel beams which 
could then be sampled by a detector scanning the ob- 
ject. Of course, in an alternative embodiment it would 
be possible to translate the object in the appropriate 
transverse directions while holding the focused trans- 
mitter and receiver in substantially fixed positions. 

Brief Description of the Drawing 

A more complete understanding of the invention 
may be obtained by reading the following description of 
specific illustrative embodiments of the invention in con- 
junction with the appended drawing in which: 

FIG. 1 shows a simplified block diagram of an illus- 
trative terahertz imaging system in accordance with 
the principles of the present invention; 
FIGs. 2 and 3 show comparisons between input te- 
rahertz waveforms and the output waveform after 
propagating through a known material; 
FIGs. 4 through 6 show illustrative embodiments for 
insuring a desired amount of scanning for the object 
to be scanned by the system of FIG. 1; 
FIG. 7 shows a portion of an illustrative terahertz 
focal plane array useful in the embodiment in FIG. 
6; and 

FIG. 8 shows an image of a semiconductor dual-in- 
line packaged chip produced by the illustrated ter- 
ahertz imaging system. 

Detailed Description 

The THz imaging system of FIG. 1 in accordance 
with the present invention includes a source 1 of repet- 
itive, femtosecond duration, optical pulses, an optically 
gated transmitter 2 of THz transients having a broad 
spectral bandwidth, imaging optics 3 comprising lenses 
and/or mirrors, an object 4 to be investigated, a time- 
gated detector or detector array 5, a scanning delay 6 
capable of changing the delay between the femtosec- 
ond gating pulses on the transmitter and detector(s) at 
a rate of a few Hz to hundreds of Hz for the purpose of 
temporally heterodyning the THz-frequency transients 
down into the acoustic (Hz) range so that they can be 
processed by electronic techniques, a digital signal 
processing unit 7 including a digital signal processor and 
an A/D converter to process the time-domain data and 
extract the desired information, and a display 8 to view 
the image. 

Certain material and objects can be characterized 
by a frequency-dependent absorption, dispersion, and 
reflection of terahertz transients in signals which pass 
through the material or object. The present terahertz im- 
aging system analyses that frequency dependence in 
the time-domain by collecting that transmitted signal 
propagating through the object and then processing the 
information contained in those signals for every point or 
"pixel" on that object. This is a non-invasive imaging 
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technique that is capable of differentiating between dif- 
ferent materials, chemical compositions, or environ- 
ments. This technique has applications not solely limited 
to biomedical imaging of tissue, "safe X-rays", chemical 
reaction analysis, environmental and pollution control, 
process control, materials inspection, fault detection, 
non-contact mapping of the doping in semiconductor 
wafers, profiling of doping and defects in laser crystals, 
and packaging inspection. 

A typical terahertz transmitter emits a single cycle 
of electromagnetic radiation centered at 1 THz after be- 
ing illuminated by a 1 00 fs laser pulse from either a mod- 
elocked dye laser operating around 620 nm or a mode- 
locked TkSapphire or Cr:LiSAF laser operating around 
800 nm. Because of the short duration of the THz-tran- 
sient, the spectrum is broadband, typically extending 
from less than 100 GHz to several THz. 

No electronic circuit is capable of measuring and 
processing THz bandwidth electrical signals directly at 
this time. Sampling techniques based on the repetitive 
nature (typically ~1 00 MHz repetition rate) of the optical 
and THz pulses can be used to measure the THz wave- 
forms provided that the sampling window is shorter than 
any THz transient to be measured. Typical photocon- 
ducting sampling gates have sampling times shorter 
than 0.5 ps and are thus able to measure frequency tran- 
sients in excess of 2 THz. No fast electronics is needed 
in the sampling technique, and only the average photo- 
current in the dipole antenna is measured. Similar to a 
sampling scope, the delay between the THz waveform 
and the detector gating pulse is scanned slowly at a rate 
of about 1 0-1 00 Hz. Thus, each sampling pulse samples 
the THz pulse at a somewhat different time, until the en- 
tire THz waveform has been reconstructed from the 
samples. This leads to a "temporal down conversion" of 
the THz waveform into the kHz range, where it can read- 
ily be processed by electronics. This sampling tech- 
nique is also known as Equivalent-Time-Sampling but 
is otherwise used in any commercial digital sampling os- 
cilloscope. This isochronous sampling technique has 
been described for picosecond optical sampling by K. 
Weingarten etai. in IEEE J. of Quantum Electronics, Vol. 
24, No. 2, pp. 198-220 (1988). 

Many, if not most, compounds show very strong fre- 
quency-dependent absorption or reflection within the 
frequency range covered by these THz transients. Also, 
molecules and chemical compounds, at least in the gas 
phase, but also ions in certain crystals, have strong and 
sharp absorption lines in the THz spectral regions. The 
absorption lines are characteristic of the material under 
study such as a water molecule and its environment and 
can serve as a "fingerprint" of the molecule. Each chem- 
ical substance hence leads to a characteristic THz 
waveform that identifies the chemical composition and 
environment of the sample. There are also materials 
that are completely opaque to THz radiation such as 
metals and other materials with high electrical conduc- 
tivity. 



In the present THz imaging system, the spectra de- 
scribed above need not be computed or directly meas- 
ured. Instead, the relevant information can be extracted 
right from the time-domain data, in a manner similar to 

5 speech recognition and processing. In FIGs. 2 and 3 are 
the input THz waveforms (dashed) and the waveforms 
after propagation through the doped silicon sample 
(FIG. 2) and water vapor (FIG. 3). 

The digital signal processor can recognize the char- 

10 acteristic shapes of the transmitted THz waveforms 
(specific shape change and attenuation in one case for 
silicon and ringing with characteristic frequency in the 
other case for water vapor), to determine the particular 
material at the spot illuminated by the THz beam. This 

15 requires training (or loading) the DSP with these specific 
waveforms in advance. Such a procedure is well within 
the knowledge of persons skilled in the art and will not 
be repeated here. 

In a particular embodiment shown in FIG. 4 for the 

20 transmitter, receiver, and optics of FIG. 1 , the THz beam 
emerging from the transmitter is focused to a diffraction- 
limited spot of 0.3-0.5 mm diameter. This is the diffrac- 
tion-limited spot size for 1 THz radiation and close to the 
best spatial resolution possible with this technique. This 

25 spot is then imaged onto a single THz detector. The 
sample is placed in the focal plane of the THz beam and 
scanned in x and y in a zigzag pattern using two orthog- 
onal, motor driven translation stages (shown pictorially 
by the x and y arrows). 

30 The delay between transmitter and detector gating 
pulses is continuously scanned by a 1 0 Hz scanning de- 
lay line. The amplitude of the scanning delay line can 
be adjusted and determines the time window of data ac- 
quisition - a 1 mm amplitude corresponds to a 6.7 ps 

35 time window. The average photocurrent induced in the 
photoconducting dipole detector is measured with a cur- 
re nt-to-voltage converter and then fed to a A/D-convert- 
er and DSP processor card. We use an A/D-converter 
capable of a 50 kHz conversion rate, and a DSP proc- 

40 essor that can Fourier-transform the waveforms at a rate 
of 1 00 FFTs each second. Thus, with this system we can 
easily obtain the FFT spectrum of each THz waveform 
synchronously with the 10 Hz scan rate. 

In an example from practice, the FFT spectrum is 

45 represented on the display screen as a colored dot with 
the frequency components of the THz spectrum repre- 
sented by the frequencies of the visible (rainbow) spec- 
trum. That is, the terahertz spectrum is mapped onto the 
visible spectrum and only those frequency components 

50 which propagate through the object under study can 
contribute to the displayed color. 

Since DSPs are used in this system, it is also pos- 
sible to utilize time domain techniques by computing the 
convolution (correlation) between the received terahertz 

55 signals and stored patterns which are related to partic- 
ular elements, compounds, etc. The signals which most 
closely match the received signals will identify the point 
of the object being scanned. 
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In another example from experimental practice, the 
DSP processor looks for certain absorption lines which 
are characteristic of a specific molecule and assigns a 
specific color and intensity to this absorption pattern. Af- 
ter each scan, the sample is moved by one "pixel" (pref- 
erably roughly the spot size of the THz beam on the 
sample) and the display is updated for that particular pix- 
el. With the above system, a 50x50 image can be ac- 
quired and displayed in just over 4 minutes. 

FIG. 8 shows a preliminary result for a THz image 
obtained as described above. The picture is a so-called 
"THz X-Ray" of a packaged semiconductor chip. 

In another embodiment shown in FIG. 5 for the 
transmitter, receiver, and optics of the imaging system 
in FIG. 1, the sample remains stationary and the THz 
beams are scanned across the sample. This can be 
done either by mechanically steering the THz beam with 
mirrors or by optical steering of the THz beam (in which 
case, steering of the optical beams causes a steering of 
the THz beams). 

In the embodiment shown in FIG. 6 for the transmit- 
ter, receiver, and optics of the imaging system in FIG. 1 , 
THz waveforms for the entire sample are acquired si- 
multaneously by using a focal-plane THz detector array 
as shown in FIG. 7. Here, the entire sample is flood- 
illuminated by a THz beam, and the sample is imaged 
onto the focal-plane detector array using a lens system. 
The flood-illumination causes the illumination to appear 
as parallel beams from individual point sources. 

The focal plane THz detector array consists of a 
two-dimensional array of THz dipoie antennas (in this 
case 50 pm on each side) which are lithographically de- 
fined on a low temperature (LT)-GaAs or radiation-dam- 
aged Silicon-on-Sapphire (SOS) chip so that the gating 
time is subpicosecond. MSM photoconductive switches 
using interdigitated finger contacts are defined between 
the antenna chip. The size of the interdigitated photo- 
conductive MSM switch is roughly 10 pm square. Each 
of the antenna/MSM elements constitutes a THz image 
pixel. The MSM photoconductive switches are gated by 
a short optical pulse derived from a beam that covers 
the entire area of the chip and is focused onto the MSM 
detectors using a microlens array. The microiens array 
and the gate pulse can either come from the same side 
as the THz radiation (with a beam splitter), or from op- 
posite sides (in this case the THz beam travels through 
the chip substrates before it is detected by the anten- 
nas). Only 1 pJ of readout energy is required for each 
MSM gate, so that a 10 nJ optical pulse can gate a 
100x100 focal plane array. The antenna chip is solder 
bump-bonded to another chip underneath with one con- 
tact on each antenna pad that carries the detected 
photocurrent off the chip and to the DSP processor. 
Preferably, the underlying chip contacted by the solder 
bumps is a CCD array, so that all pixels can be read out 
sequentially like a video camera. The photogenerated 
charges are accumulated in the CCD array over many 
optical pulses before the charge is read out. 



Claims 

1. A method of inspecting a medium having one or 
more absorptive, dispersive, or reflective constitu- 

s ents, comprising the steps of: 

directing a sequence of high frequency, elec- 
tromagnetic pulse signal at said medium, said 
pulse signals having a waveform with a spectral 

10 bandwidth sufficiently broad to encompass fre- 

quency dependent variations in the absorption, 
dispersion, or reflectivity of at least one constit- 
uent of the medium; and 
detecting the directed pulse signals following 

15 one of reflection by and propagation through 

the medium, at least some of the detected sig- 
nals being characterized by temporal distor- 
tions representative of said frequency depend- 
ent variations, to thereby obtain spectral infor- 

20 mation about said constituents in real time. 

2. The method according to claim 1 , further including 
a step of: 

analyzing detected waveforms having tempo- 
25 ral distortions in the time domain to obtain said 
spectral information. 

3. The method according to claim 2, further including 
a step of: 

30 creating an image of said object based on 

spectral information obtained during said analyzing 
step. 

4. The method according to claim 2, further including 
35 a step of: 

identifying at least one constituent of the me- 
dium based on spectral information obtained during 
said analyzing step. 

40 5. The method according to claim 4, wherein said 
identifying step includes comparing spectral infor- 
mation obtained during said analyzing step wit 
stored spectral information associated with a plural- 
ity of compositions. 

45 

6. The method according o claim 1, wherein said di- 
recting step includes causing said transmitted sig- 
nals to propagate along a plurality of spatially dis- 
tinct paths through said medium. 

50 

7. The method according to claim 1, wherein said di- 
recting step includes causing said transmitted sig- 
nals to propagate along a plurality of spatially dis- 
tinct paths so as to be incident on a plurality of spa- 

55 tially distinct areas of the medium. 

8. The method according to claim 1, wherein the trans- 
mitted electromagnetic pulses signals are in the f re- 
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quency range between from 100 Ghz to 20 Thz. 

9. An apparatus for inspecting a medium having one 
or more absorptive, dispersive, or reflective constit- 
uents, comprising: 5 

a transmitter for directing a sequence of high 
frequency electomagnetic pulse signals at said 
medium, said pulse signals having a waveform 
with a spectral bandwidth sufficiently broad to 10 
encompass frequency dependent variations in 
the absorption, dispersion, or reflectivity of at 
least one constituent of the medium; and 
a detector for detecting the directed pulse sig- 
nals following one of reflection by and propa- 15 
gation through the medium, at least some of de- 
tected signal waveforms being characterized 
by temporal distortions representative of said 
frequency dependent variations medium, to 
thereby obtain spectral information about said 20 
constituents in real time. 

1 0. The apparatus according to claim 9 further including 
means for causing said transmitted signals to prop- 
agate along a plurality of spatially distinct paths 25 
through said medium. 

11. The apparatus according to claim 9, further includ- 
ing means for causing said transmitted signals to 
propagate along a plurality of spatially distinct paths 30 
so as to be incident on a plurality of spatially distinct 
areas of the medium. 

12. The apparatus according to claim 11 , further includ- 
ing analyzing means coupled to the detector for an- 35 
alyzing detected waveforms having temporal distor- 
tions to thereby obtain said spectral information. 

13. The apparatus according to claim 12, further includ- 
ing means for storing spectral information relating 40 
to a plurality of compositions and wherein said an- 
alyzing means comprises means for comparing ob- 
tained spectral information with said stored spectral 
information to thereby identify a constituent of the 
medium. 45 

14. The apparatus according to claim 9, wherein said 
transmitter is operable to transmit pulses in the fre- 
quency range between from 1 00 Ghz to 20 Thz. 

50 
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X'h 0 . ZtllZjZ 10,^11 S-ftA-TS . ?^2r# 
[»13] 

4 \ 2d ) 2 



I=N 0 e 



■<t 



[0045] i^CiOs 5-f 6waj*«SE 
1 2 f* ft . m^«S03t^±#' *) 1MB d/w + i x« 

{4. ifcsreasfu KpWJl, mffiP^sg®ido2 

S^JtWTfto 
[^14] 

^ 3 



!( 7) 000-135268 (P2000-135268A) 



[0046] ^tihcomimmi 3 t&a-t. ®v^*jet 

KffiM^<tS«±, *IP H lgBtid£lmm. 
J± P £ 1 0 a t m . EPJn«ffiV^ 1 5 0 0 V^fS^-f- 
4. £;hX49. A^itfOBSISItSrl OOjusttS 

[0047] ¥#f*5^m^a 1 6 b(4. H 1 4 HZjprtX 

■&IBfc: X££#Alt^£ „ *awx£f&£-r «■ . P - 
n$H=r<9t2j&\ P- i -n^:MT<9&c0i>WrS 

c: & . <i 4 o ^mimtamm 1 6 b zmm 
[0048] -mm^-i yty^i 6 tzx 

;1/*K SiaWfcJJB^^fflfcj&sfiqtfS. £^4 3 

0 &stEH^PB£S£-rs(~(4. — «wca. -e^v 
•ess. s^^'fe. ^rawt^-^aufi^) 

[0 04 9] y-fe>"»f 1 6(4. 01 5t^tJ:3 
figgl2a#4ft3i (x, y, z) (Cfcl^K 

xttS-ftKTfPSra-a " ithC:@ES£*l-2>. 1 2 

a rt-C*0Wj&*£3r 9lir~>O^I>k, «K0<Offl»^tt^) 
WJ#fflltT'&£rt;#>. 01 6(^ti^:, *#tfuO 
filC, £$ggl 2 afl<^*J£4M^4#*^H3 2 
SrlftftSifc^SStV^ d<0^S^S32(4s -ag© 

1 2 ac7)T^-WIS§nfc2 0C7)X^>-i ;: n-7 3 3 
a , 3 3 b k . 1 2 a (TXtfrfcEBSfUfc 1 

«yyn-7 3 3 c k^fl§j£§ix.. Zti^comizQ 
M12a fca^-CKStcSBiS-frS £ k *C*flr-3 izM 

32(4. 3o^yy'n-7 3 3a, 33b, 3 3c 
(dR£-f, *ftja±<o»o*j#y^n-5T«J« LX h 
4^„ 

[00 50] =jA V-fey-f 1 6(4, ftHgl 2 a« 
Ztl&lZ-otlX. m&MSkl 2a<7) P 1 , p2, p3, 
-, pn^xaURStsK'f^KDajiJt^fBi-S. » 
^yhp 1 , P 2 , p3, — , p nc0f>yf-i4. * 
itt^l^, 017 (C^cf 4 3 fc. 0 . 7 6mmt 

[00 5 1] 01 7f>J;3fc:5>f ^-fe^m 6<5D*ftbK 



-fyht* 7*;«1 l/o > & *£\ ^-f y-fey-f 

7 7-C5S£*U>. 

[0052] 01 5^4 3 *l7*)VMb$m 
jfiftfiSSflTV^faSBRl 2 acoaiw^ y hp 1 , P 
2, p3, p4, P 5OX^3§*O^7/Sr01 9*» 
£>2 3fcijct. Srfc. 0£fgB&C/t&^#>. 
(42^ybfe^(^$tlTV^. 02 2*>4>2 3(Ci3 
Wt, W>bP 1- 1-C(4, ftiSS 1 2aco«:J;5 
XH^3i*««s^ffi$tl. Wybpl - 2f(4. 
tfT-feJUOSBSStci £ XIS«Ba*«0SSfi8B## { Ja**& 
A. Wyhpl-3Hp3-2 £X*teAy°-bWP 
AWtCfc******^^* . 4fc, W^fp3 

- 3 j&»6«tf@Wfc: J: -S XttSg^^SKStfi-tfJL^W: 1 
«>. ^yfp4-i -cti^r-fe/nc J: &x«3SBHi<o 

««#vh£ < 3r£ . y h p 4 - 2 frhOtJ y h p 5 

- 2 * xnmvuz 4 h xnim3:comm<vfr k & o . # 

A y h p 5 - 3 T14-a^S 1 2 a 0»(c 4 £ X^3®* 

^ff-^^3i^7C^7 7^-fk. 024<5D4-5(c^ 
S, 

[0053]^, MK02 5 £^lT4 3 tS^S 1 2 a 
V^»^i. fflltU, Xifta3Mt^tt[£BL. ^coS 
M^Sk. 02 6 (a) Ki2 7 (d) 

j&^feii*. zcomizts^x . *i«00 0 0 0-1 1 1 
K4. 5^ytyfi 6®iHiiSfiit. ^W(4. « 
aas-cs>i>. 

[0054] mz, xmm&mc?)mm?-?frt>&mgk 

1 2afc4l>««»^M«t. S^ty>fXK^x v 

rtj»x«t5Hi<oi*auiL*>«mT'S*v^'c, m&d 

[ 0 0 5 5 ] 02 4 v fWmt*. IffflOjfO 

XKSUDW &m*«(c 1/2('LT«J1LTV^#\ 
(CJiJH*^ X«»J|XSS*n»|45 tffliW & «> y W77[6] 

tzmz, mfrm£®vxmiEtz> - 1 & , 

[0056] yWft^li, «W«0¥ffl!l*fC*4< 

%&izt*<Dm±%:mmxn^\ ^commit, tsmzti 
tzmmtmxhhMm-^m^r—^^zx otwas £ 
t ^'t- s s co-o -as § *ife-c* * a mmco x«ir* 

X XHS5iSO«iaj-r- ^ Srlil LxmKf-? *fp* 
-f'-^'— Stf"-I> ^S^Sr¥'J^-rtL(4' 4 1 v&» & S . 

£iii?>**-R l*^ k v-» o £ fc (4. ism vtzmwmm 
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[0057] ^Wtl, mfeWg?-? k LT02 9^> 
H 3 0 «t 3 ^HSIt"- * * -BrT tJ S , 

Igfij^a^x-^ kSffcaj^^&^'-y k **t« LT 
fJS-T Sti?**. 3&7C7 i '-:?k tTlBISt^t^'W 

2 9 {ijtL^c7)^M-C04 7 tc^-Txy h -;K H3 0(i 

H^itat r a j comtfmwztLfzmqx'm 4 7t 

[0 0 58] ?-^#W17(l 03 liZ^Xd 

*)V¥3 5^LTS»W3 6lz3m-f&t fcfcfc, y 
V - A 3 7 IZM K) Wf h tlfz V 1/ y -f K 3 8 V 3 
4 ^SgP^S^ 1 2 a tcff Uf+ft £ AaCLfcfeW 
SrffifflTJ-S .Ii:#'T"£.£.. #-^y34i4. V W>f H3 
8^'K)f^t^:Vi>B#^. -S^l 2a#>£>Mix&J:d £0 

usnu o fc„ y u-j* 3 7 tjg(estL3tai3 9 tjuow 

[00 59] fMffllJIE 1 8(4. Hi {ZTfCfi. 0 fc, #36 

y-cffi^-r s i k * s t-# zz x\ mnmm. 1 s £ 3 
2^g^iWKa^«7 ; '-^*iMfi-rSo 2^g 

[0060] *ij»§gg 1 8 te J: SaM^fcitX-r- * 
T^Hifi^«t=*3 v •if , StfiBH 1 k M»m tztt 

j&tz mtiz\$m—w%*tt lx . wR&wr*. 
[0061] <m2mmw>m3 2(4. 

irctt .010 <d mriam 1 mmmmm\m&mu t n 
Mhz, x.wmi4ty4y^y^riemm^fix\^ 
#\ ztiizMz-x&mski 2 a ^^§n>taK«is y« 



Jj&]frt>mB-f h CCD^740 1 2 a C0± 

[00 62] ;«ccd^74 OfcJ;^. mqco&M 
r - * k HiWtfeT - * £ JDtfW" S-k ^'T- 1^. cc 

WWX K) lEMfc^r h . XUft-r - ^ "Ci40 3 3 t^rfi 
fflBC CD^540 TTOHJ Lfe*HRTR£f* — ^ Srftffl 

k * < ^facomwtfmftx* § s . 

[0063] <lf§3^y&M>ll3 4 {4. *56^<OH3 

5- C14. y«*l*ltilSS*ufeXjR«l 4 k 9-f >-^^ 
-9-1 6lztax.X. xW^-[6)tt>X^14k7-fy-feV 

1 6#lEa§tlTfcD. fiHgl 2a£9 0° -fifut 

2 ^ri*]*^ xH*»RH*fc!flBt3 *x S i "5 fc* -> X v ^ . 

W^-iejcoxiS^Mo^-^^ior. a^r^Tfi^^ 
tt^iSifiJ^TOWS©^Brtgk=5rS. 2^[fti<0 

3 &7cM££fl^£* £ ; k ifiV S S . 
[0064] <^4H»E2i>ll3 7 {4. *%BJ<0m4 

{4. «B0E#cojaHKaH£fliil4 ld«KStiTM^ 
ttSS4 1 ±fc:x«« 1 4*^«i™fc:««Sii4 k k h 

awfitt.514 1 ^hffliJi H^iS« 6^" 

ESStLTV^. S*c«lil4 1 kSStt^y-fey-9-l 6 

fljiBfcyi, 2O60 u yyvmuz* u -y h 2 1 m^xts: 

6- ?X7 2 6mm%fLX^&« Xlf 14i»1^4ft 
OfeS^l 2a<D|BIfctt. T)V$.—*72*frt>-%hyi)V 
■?2 5 am t hit. XWg 14k —WlZft-o TMft 

[0065] mmmimmwizt5^x . 03 7*^0 

4 0(;^tJ;^:, tailSl 2 aSrzW^^iHSkk 
fcfc. X^«14^FJTSftJe^Sfe0S-«i:-g>-k^4 

2 a o»as*£«»: l , y a x<m& . x. >y s^mm 

04 1 fc^^Bar-^^feix* o 
[0066] <lg5liil>l4 2«. *%B^|p5 

«4. m3Abm3i<Dm^^t>^->tzm^^Lx 

V^S. -T5r*>t>. 9 0' -ftif^r&]^2^0X$l«14 
^UaSiXS k k fct. mft&faX'-fo 1 2 aco 
JHHt«rlB2o»X««l 4«ii«ii3l4 1 

(4. 2-3«X3f«14i490° COfenfflffi* 
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% h 1 t i>iz. msmk<wi-®mjt*&i}tofrt>mfc? 

[oo67] mu^mmmmx'H. xwm 1 4« i o 

• t y i-commm y h xxm^mmv^^ 3 6 o • 

EWRU -S^l 2a{ifriBl 0° cO*lt7K-f 
SStO. 7 6mmt"7 f-T' z tt^ltfj^^ifrf ft i fc # 

[0068] <memtmm>m4 *hbjco^6 
TJim 4 2 § tv/nsfflx&> tjux-mm vx v ^ . 

AM7/H4 2ii > R*T*ffiffl§*lftfctf)-C\ 
ittl>< : 0?:R5lt-rS.5t^cO7°vXf->y^^^>yr4 3 

•4 7VWK4 2J4. 04 4 KiRfJ; 0 fc. JK5£^l'h4 4 

mzm%^tzimxh%> o t . mummmm^m^ 1 

[0069] <WStW>&^ 03 2 fc^-TSS2^fi 
H4 5, H4 607n-f-r-b^^TiM^-ri»* i \ 

wmimmz^ wc *> t star nmx'$> ft . 

[0070] , SJjes^fcovvrSHlH-S fc . 04 
5£*-tJ:o^ Xf'yTl Ol-CSSfil 2a£>zW 
*I^I<0«SUi«* U -fe 7 h t . Xf y 7 1 0 2 
IHftW*. Xf771 0 3X j mgk \ 2 a^ma^^ij 

jtusavfcWBrU mm-tui z.f->yy°i 0 4xxm^ 

* £IB«-f ft .»KsXf77*10 6 T- 1 fS#tf>«H 
<0£T <0«aj jftTOr - * JR»*^T t it * mn 

TW-f - ?K##'5£T LT VtfUf , XfvTl 0 7 "C 1 
ta#<9*KH&CCD#;*7 4 0-C4RKU Xf-y710 

fv71 0 1 l/ZM 9 . SrfWfRTf * . 
[007 1] <f-?)M^^, mff5«!l5£»ff 

izx ^m^Ltz'f-fcommm^mA ecoyv-^+r 
-hi,zv£^xmw~t&» 

[0072] Xf>y71 1 1 X\ lEHt/t^^ > 

-fe>"+M 6<affi:rt7W*R*tfJ-$\ Xf 771 1 2 

f 771 1 3T\ iO^ifSO^t^-feV-fMSi: z« 



SWcfcHESifTS^JS-ttS. Xxy71 14T14. ^ 

tbim tfz xwtmmzm^ , 5 4 1 6 

jw=, xf77i 1 5x\ mmvtnMi&M&WBW- 
ft. z.<D^>mmi±. mcDXT-vyx'&MT-fbMfc 

[0 07 3] Xf»/71 1 7T14. CCD*^740t 

1 1 8t\ CCD^54 0(7)feif^};7^y-fey 
-9- 1 6c7>X*iaji*7 ; '-^ Sr^PJcf'^MB^^-tirTX 
U%mm-&, Xf'/Tl 1 9X14, HUlBXT--y7°l 1 
8 VfftR V fz-mtT- ?frt>3 &jtW&r— 9 * ffifcfr 
ft. 

[ 0 0 7 4 ] Sjiivt, Xf7712 0-Cti, WS'J**X 
-^SrK^ffi-to ^TOiMiJK^JS-rftfJS'J 
^flt^ai LT!H^MS&t ft i 0 *> , 

LTv^ft-a^®i 2at-ans<is^&«sfUJc*fie^-s 

W&m*lW£&< 'fro ZbifiX't &<VXft£ U\ 
Xf771 2 1T14. 3 &7CiJ#-f - ^ «»^J4"^fi[S 

«^^^7>-h Lxmrnmico^mm^L, xf77 

1 2 2X14. MfSXx-yT 1 2 0 "Ctt^ajLfcfWOBIS 

f 1 -^ tttmz^xmmnwmzmML. xj-vti 2 

[007 5] Xf7712 4T'(i, S^12aW 

^ft tzMz . WJ»m-f'- ^ liZjttJS t^rv ^ftm-r - ^ 
*&jW5*»*WBrt, &tt{4Xx 7 7 1 2 7 tWn^ 
-?^UXf771 2 8tWfl^*l> 
Wi'Xf771 2 5tCjStf. Xf 771 2 5 "CUt, 

VW^t^-iSc-rft^A^^K L . ^-Sc t=5rim 
Xf 7 712 7 T*^^.ftV-i7 Srft-^-LT Xf 7 7 1 

2 8 T^iln a n^KSr t . tT V ^ifXf 7 7 1 
2 6f^|g^t6. m±WWSIS*{4. IIMff 
-^-i:MJSB§*Tf2'[S^ft. 

[0076] <^a^x-^^S>^.Sii: LTWJS'J 
§tLft«BW4. OT<Ofc^* s fcft. 

i . -a^^^KAfiaofra^xva^ 

2. 

3. g©w<7)m 

4. ^fOO^M 

[0077] nnB^Jft^Kfrttiti 
coi,h^#»]L, %&^<t-5T{4^^^0it-r^^ 
&ft„ ^n a 0 cr)#»fc{i. ^^^5^7, 7if- 
^>*L**», -aSS^O fOSP^t-^cOJ; d 2rT 



(110) 100-1 35268 (P2000-135268A) 



g izttm utv- ? & & <r k * u ^„ 
[0078] mt?MM^&<7> o &•? . -ass^wtx 

wkwm-h. mm. hhm^z^xmiAiB3 

HftAf^g- 1 §ek i^JBSB^lB 2 0ttAft# 1 #5^4 
rfij , r ? j milzW^Zii. r ^J tm^titzn 
r^j t ®¥Ztifz&MMzlZ2m*-&M-£ti-%:lW$ 

mm&t&b. r fflj t®¥-£ti?z&%iMt,zi±im. 

rfij k £\i¥$titz&Mmzl±2mb% D . EP^$ft£ 

[0079] mm. ze>*o ^m^^tt t . 

iz. -ei0 2@O^&a^ilS^^^A7]LT-a^ 

1&fB^m£¥'MrSk, rigjj t r^j <02S{:S«t 
£ Mfcfr-r- 9 &&%mWV>M 0 &##»ffi ffi k LT f&S 

[ o o 8 o ] yac pft^aA, m^wnxuxx/mm 
<n*?. m^m&^-fizm&L^'bmMbte, ±xm 

%SCOm\b LXmttZ . 062 (A) (B) 

\,Z^f ± 3 fc. ^ffif-^* s IIf-^^9 0^°-- fey 
hk-ft-fStf, 1 0%m^^I^'#^tTS-^L 
£vtffr&. &ShbLXm®LX. ZtiZ&*.X*H Xtf 

4^^-^ tmmLtz^n. mwmtit txmm 

WiZ. H6 2 (C) \Z^tX.o\Z. W&f—?tf 

mm?-? iztt Lxm&miz-wcth <,zi>mbt>-? . if 

-f Xj&*A* ^XI2>h£ v^>«Sfe&»fef-^«^— *4: 

ttJsELTlE'SLTfc*. Wr"-***, SBttL*r*3& 
HoKHfc fctitf , Hfcfrfcsa L JtagfflTJb & k WJ L , 

ximtz. 

[ o o 8 1 ] tfz. mMmmmbwmmm<?>im<n7F 

BkjmnLxmz&ar&izn. $m#*v hcomms.* 
tfmmx°fo&b3§r£t>ti&cr>x*. ®m%&±lx^— 



[0082] <«Wfttl»f f --*>«VvC . mfflffl 

-?it.<!zco3 ^mmxtm Lxa<z.bt>mtL 

3 . 4MKW^9-Wf8**tf 3J>C7CT-^1f-f Xi: L 

m&mtthc\bi£Xlbh> 
[ 0 0 8 3 ] 04 7J4. Mm?— ^7CX- ^ ^ 

MLximttzi>cDX'h*) . mi o^mmzn^~thi> 

cob. 03 7l l znifc-t&i>cDbti i ^tiX^Z>« 

[ o o 84 ] si o cr>mmznm~fh ? « . 

{±02 9, 03OOffMT"fB1t£;ft.. C^IS'lit/Sx- 
feA. 01 0<7)§ggfc*fJ(frr& -r-^ti. «»tXHW 

mure* =5rv»3&s, sepjtttsaiaWHK^^rtisa! 

[0085] s/s, 03 yoiisfc^je-r&x-^f*. 
S fctf>, Kffi« Jt=5r if^JSiJtJ «t ifftm^mm: b<?m 

m<o*-b vxmm-t&cibtf?%&. imt 037 
Kfiwi hfrbx-h 0 s 04 2^04 3«oiiirrfc« 

m-fhc\bifiX*%h. 

[0086] <m^com&>mmmmmtmkc7)?-5> 
a. mbLx^-fizmmzti. mvto&mmmzttm 
lx. &mmTLtzmi\<owMm%i-&c\bx\ &m 

[0087] c\<7)b%. xmFTX^y LfZT-f 

it. xmrnf-fx-fo*). iox-^Srit^s^ 

[0088] mmmmmm^-fiz^tbm 
-ft L^^imf- ^ s- , mm&b-wct&sm-T 

-mm*. mmmzttzzfomzmiLizT-'-fiztt 
-5- u t>mm&>mwmwx'% z>£o iz-tz, . 

[0089] aftH3WKO|««±, l*]ii«jtSrifc 
^-f 4k Ware* Sfeft, iM'JDfffi«««3tSrmria^ 

k -iwt^^-rs z\btmtt^\ z <m& . mm^ 
mmzm^thtz»). MiPix?-*?. mimzwimmxm 
■&i.xm^-th. #5>-m&-im. yixwftfoifi 



(?l 1 ) 100-1 35268 (P2000-135268A) 



[0090] <^J3!|3- H^)BIK> 3&7GM»fc:lf£$ 
«03-K&t)Si8W-«.ii:* I 'C*-6. Wit** H3 0^ 

[009 1] CCD^?^7 Srfiix. Jt^STIS . ^ 

c k a s a6«afe?^ES'J* s # < , c c d # 

* ^TNOJRftT*- ^{±-?-c7)fe^»Kffi*t & £ k "C* 
-?b UOR»9&tJ*£\ *V-f MM»Sr^3*\ IK# 
if s C CD*^ ^T-^fLfcT-'-^WJx-* k LT 

a, # * 5 e \- y m%WM. ? ■< how h § 
2 . s*afe& i mxmrn u » i aiii±«»a-fe 

3 . Xf772 T'leit Ufc 1 8J2Lhf3*3S6fe'r - ? * 

liitu^ts'ii l i sjii±^**a-fex- * t am l . x 

T7/ 2T-|a«L^a2pafe7 ; -^*«iEL^ 0. 

4. ^sjs^Hff|gXT--yT2. 3Sr^03i-r. m 

sett t £x-^o©tag#'r&]±^& . 
[0092] £*l<->tf>Wt\ mfflco-&MT~?-*?mtt 

a-v commix # *u* . <r ^flHB^Tcfciaspjcoasffl 

[0093] minmftfr- H«BW4, OCRTWI 
[0094] 3&7tWH\ C CDW ycDUmZ, A 

[0095] Z\cr>Wr£. mSHZ^At LX . 3&tcH 

b LXM?>-M&Zb&X%%:^> 



[00 96] tfz, CCD^^7»ltf-^^t 

[0097] 3mmw»n&\±. j&t Lxmnxz 
A y*$f~?x^&, 

[0098] z<7)?j vmmmt^ ^ i/- v<mx° 

[009 9] Mi5>f v^i^S'Jt. i<^> 

b t> tclSfi-tS i k* { Mt U\ 
[ 0 1 0 0 ] a^HSr^^v^^ffifgcT)*^-, 

jtroc thlmj buwrmm^zm^tt. mm 

[ 0 1 0 1 ] i» X a 55r^M^^-Ct> . 1 
aBPfcS^EMEfirrtftWJW"* - k *fWC*4 . 
^Sfc: «t fitf . CCD^5 OSMR'r - ^ co^T- «> » 

S V \h C0X\ m#^r O C R V 7 h TfiBlTS s $ 
tc, fiJS'JCWSkWBl^SL^v^ 
[0102] ^it. -X^yfA X^ s * > b*« 

[0103] OT. 3 ^jcH^-r — ^ 60#JS: . 7 -r 
[0104] <3<JC7cW«7*-^<0ffjR>H3 4^04 

jjT-fZWfcftfato ■ziffiKv&'ommzmmz-txffi. 

[0105] 1 0 ^)Xit"-A«i[ii* (4BHI) cort 
g|$£iffilrf &k, 4SW<0#JfOifi^. 

[0106] ico#Jl36»fe»j£SfLS«0B<?5ffl(£**3- 
--y htt^^LT#X, *iX-? - tl. ul, u2, u 

3, u nkv^jfisjfiaac&fcorv^k^k. ii 

[0107] 
[*15 3 

/ = 7 0 •exp(-« l A)-exp(-u 2 A) --exp( «„a) 
= /„ • exp{- («! + m , + — + «„ )a} 

A«i#fi[J---y K^StftS. i<05C*S9g 
[0 108] 
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im 6] 



[o i 09] friai 6^«i. UiA*#3.-7 h^as 

[0110] H4 8tC***-J: dtC, 4ffltf0^3&^^S 
ffl«WX»iKKffi|R*-€-^ii. X lt X 2 , X 3 , X 4 fc 

u ^ti^tiMcom^yii = i, x 2 =3, x 3 =2, x 4 

[0 111] 



A 


Xi + X 2 


= 4 


B 


X3+X4 


= 6 


C 


X!+X 8 


= 3 


D 


X2+X4 


= 7 


17 


x x + x 4 


= 5 


F 


Xg "h Xa 


= 5 



[0 112] U&»U d<0^rffi{±H^fflVM>i:. 

[0113] WiEIBjas^#i» LT . :SH5Ji<9H3:-£ 
[0114] 4-f . ^#:OT±^fl^^S»^*^« 

49(a) iz^-tXdiz. ABj^fa<7>r-*£tt;«L 
t:b%<7M£M. -1, l*2^cr>'mzmftLX. Ztl 

McO^sE* k s H4 9 ( b ) tH4 9 ( c ) Oj; 

F{p,p)- JJ/(^>>)exp{- j2ttp{xcosp + ysm ppxdy 

-- JJJ/(jc,y)y(xcos p + ysin P - R)exp(- j2xpR)dxdydR 



[0115] fli2W»i£i: LT. 7-'Jx£3ft8#*S. 
f (x, y) «7-'JX»i, ife^T-S^tU.. 
[0116] 
[»17] 

v) = JJ f(x, y)exp {- flMjx + vyffixdy 

[0 117] ju = 0(0b^ii. &.&rC I 3-X.$>1xh. 

[0 118] 

HR18] 

F(0,v)= J"(f f{x,y)dx)oxp(-J2nvy)dy 

[0 119] ZZX\ g(y)=Jf(x, y) dxi 
nt v i^g (y) (4f (x, y) Syi^t: (ft 
») fltfl-LfckOTfc 1 ^ 3S18l«-<01<K7C7-yx 

[0120] y «rtrift|0>K#f*-* Sr 1 

7-'JX^iUct^*F (0, u) k%-oX^&« Z. 

coztfrt>. ym-fifacomZT- 9 IZM LX t> ffotf 
[ 0 1 2 1 ] H4 2 §SStf)*£\ m5 0£tk 

[0 122] 

mi 9] 

G„ (r) = JJ /(x, (* cos <? + ysin 6> - R )dxdy 

[0123] ZZX\ S ( • ) fct-r/i^raft-c* 0 » u 

11: xcostf+ys i n 6 l = Ri44 1 ^ **$><9 
IBtt^RT* 4»IWSri*!I fc R L . g 0 

(R) (4^<!^JLh<^x-^*W^t/S«>0*Scr 

[0124] ttz, 7-yx^«si 7f4. mmm 

[0125] 
[^C203 



[0126] -tfthh. /3lTfacoW3& ) ± 

<o 7- yxjHfcttx, ysiatcj3tt&-?-<o^rifti<oiea<!o 

[0 127] 
[»2 1 ] 



/-(*,><) =|[F0,v)exp{;-2ff(/« + vy)}rf^v 

= £ IJf <i(9£°F(A^)exp{/2>jp(*cos(9 + ysin £3 - 

[0 128] ZZX\ M4 2lZ^-tmmz}5\^X , S^j 
mco&m c f-9CDm$.Bbe (x, y){i, 

[0129] 
[^C22] 
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b g (x, y) = J g g (R)s{x cos 6 + y sin 6 - R)dR 

[0130] mmmz. ^fth^e-nn^^n 

[0131] 
[$C23] 



[0 132] ttz. se (R-) tis e^m^.Lx^th 
n*nb^hvmxh&t>K zny-v 

X.mmTfo&F (p, 6) Z&3£MLtzi>n. tth 
[$t24] 

J— co 



X^-thfth<T)X\ 
[|fc25] 



[0133)^25 b^2 0 ^ttr^lK, S20 <?M. 

^2 0^'-?£**l>^6(;r{4. F ( p , <9> c7) 
ffcb9£* I p I F ( p , d) SrfflV\ftJ4 <fcV^i^^ 

[ o i 3 4 ] y- v^^mx-^mi-mimmrsxn^y 
y-yyv yymmzz. 9 . y-v^mmxnR= 1/2 

[0135] LfctfoT, 7-UxraR<03SSa«*«» 
5 b ( i5t), &5&&*#£>*t& . 
[0 136] 
[R2 6] 

= — sin2?sRr , 1 (l-cos2^Rr) 

2(ar) 



[0 137] r = nw. rr 

[12 7] 



= 0, ±1 , ±2, 



g(*w)= 



^w 2 



*r = 0 



Ar:fB$fc(*o) 



[0 138] £*U± % H5 1 tc^idJSrBBRk^-S. 
H5 2te*^-J:dfc, 05 1 fcBcaUsHRfcft 

[0139]<7^ /P*Bra>aBW8aE-f-*tt» WJ 
[0140] MS. |ffiX«t«ffl#«T«llBSn&f f - 



[ 0 1 4 2 ] &iz, XT-bfrKblZ. rt»* t Wa*if<0 
[0143] ^h@ISi^tt*«lfmW^I^J{4. WS'J^h 

[0144] tuiB^at i o^as t^wajx-^tt. ^ 



St 



[0141] *Ngi6asw\ 



[ o 1 4 5 ] H3 3 mnuzxwtzmMLtzbz&n 

m-&. i <7) 2 otftffctt 5 4 y {4 s 8»ffi#<0lWIJgP 
03t»S*L4^ td-^rcoBH^ ytt*WIHc3fc$> 
t*rt 2 o<7)313I85 >f y*qSBf$<lS*»i: 2ocD3fciif 

[0146] iO**r<7)§I$#i4. x**fct*ajlStf)ia£ 

[0147] a«tai«s*5-f y-fey^Htssa 
mm*-? itwsx-t hwm. ftmm^tf 1 «fcia< =&■ 

XMI^til 7 4y^ywftffimm.it±ifh (f 

[0148]«fc, Eafffl<0X«CTX^f'Vri:3|s5WH 
<, #>o s J:fctmMi:"rSWJi»JK*«<K^r«»^ 
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com * feTOiJT" § & *> OT* < T J4& io V > . 
[0149] #&Hji£:a&6&ftSi*iSli, # 

[0150] wge*fl*iWM-4fctt. M^m*^xn 

[0151] ¥»8KH4. atoffiflMOJ***. 3tji*T" 
T'#ftV^. 

[0152] B#r- * offliE^7c(i . B«j<o£J&^ 

A^jsgtss^ s u-t^^mcomm. j a x<m. 

[0153] 1WE¥i«^3«eaSl^jE* 
[0 154] «K*tIE$-S*a\ #»k¥»0*li^£¥ij 

M-r&^^&o. *<o##atfH4. tkRwowmg. 
[0155] ^^§r^< ^iir^a^*^ . * 

urns ztx\ *tffl$-b wafrwiuf-tmi-t 
zbtfT^z. micomm^mmcoM^, m$ 

[0156] ^K^cowsmmwrnt, mm^mm-t 

[ 0 1 5 7 ] El 5 3 14. ®«f»l^i|Ltf-^f-/ 
[0158] MMf— ? tf)ffi5&*0<9 4 a K#KLjfcMiR 
m.if 1 > X<73*£\ $m$kcom$: 480tt, B«f - 

^cofii«^^«wcffl^s. 

[0159] SOSiOW'Wriis K£Lfc5ftKKB&t 4 
OfrW 05 4 (a)K!5 5 (g) fcjjr«HHfc# 



BfcyBfcTfflv^ikj&^S*. 054 (a)«M« 

a^capw- * t , BflBt 4 o %m § fufcaat* * . 
ajEMa-r -s. . 0 5 4 ( t> ) <oKisi«£ffiffl-ti> t . 

fU H^attT^SlSil*. 05 4 (c) <9K#80B 
05 4 ( d ) ^»X^««iH>05 5 

(g) <mnmm*zm^z>b, mmcomm^mm 
izmmm~t$>zbtfx*&. zoom. 05 5 (e > <n 
h^mmmt. mm^mth^, 055 < f > 

[ 0 1 6 0 ] 05 6(5, Cn^OMSC^fc^jSt^AliJ 
^1P«^^ Ltzi><DX\ r W«Xl«f-^*ia 

[0161] XH^tii^a. y-fey^T^JiJ 

[0 162] 9>f y-fey^-cttajb^aj^jfc: 

14. s4X&£&ti&. z\<n>JAX\±. xmcoiRMM*? 

[0 163] zzx\ j a xzwmtz z b x\ mmco 

[0164] m-^ffi(4. JWtTWft-C**. ffi^lX 
tB^jX^ty^V •> k . J A Xifim-oX h MMWfc. 

{43 3 ymn^x'^-i t>tih <ox\ d&x*^ 

[0 165] 
[»28] 

g(')= t,f{k)h(i-k) 

[0 166] HRK:, ¥»7-f;^tt, h (x) toe 
^ s cf,^ ( x = 0 ) (DUftX'Jz* < . +*»*»4>«t* 
fc«vvjN*<»je^-4<o , C, M^fc^tSH <±W) SrS 
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[0 167] 
[*2 9] 



fl/9 1/9 1/9 1/9 1/9^ 

1/9 1/9 1/9 1/9 1/9 

1/9 1/9 1/9 1/9 1/9 

1/9 1/9 1/9 1/9 1/9 

[1/9 1/9 1/9 1/9 1/9) 



[0168] JPM^B«£^%1~S%£\ h ( i , j ) 
[0 169] 



W 



[R30] 

(115 0 
0 1/5 
1/5 0 
0 1/5 
1/5 0 



1/5 0 1/5^ 

0 1/5 0 

1/5 0 1/5 

0 1/5 0 

1/5 0 1/5 



[0170] m*-tt$ s F*nb-t&Mte&- <k<v£ o fc* 

ft. 

1/8 1/16 1/8 1/16 1/8' 

1/16 1/8 1/16 1/8 1/16 

[h]= 1/8 1/16 1/8 1/16 1/8 

1/16 1/8 1/16 1/8 1/16 

1/8 1/16 1/8 1/16 1/8 



[0171] 

[R3 1] 

0 1/6 0 1/6 0 ' 

1/6 1/3 1/6 1/3 1/6 

[h]= 0 lie 0 1/6 0 

1/6 1/3 1/6 1/3 1/6 

0 1/6 0 1/6 0 



[0 17 2] ztit>yj>u?comi,z^ 4)V9ifi 

[0173] y< frfffi^focowmomMMzm 
\,zwbxwb~?foh. t.iz. *?4 tv~7 

[0174] J A XM*^vflSO^, ff!57 4 )V 
9 *$m*h <I k T-x -y i: % < J A X 

ttmz$hz\ki$x%h« 

[0175] ZCOX 0%7 4 frfym&BMiZ. X >y v 

m^wm^ttzt^x^-mt^-t^mitth z t 

X*jm k & K) , ^fifclix >y >?gp^2r^ffl-f- ft^^ft 
[0 1 76] v^4. H5 7tfc^T, B«*fc«3l*4 

a 2 =j(x-/i)v(^=^, 2 i-^^+^^u-^y 



[0177] £<3HflW»te»W'HRWfcfc^T, 057 
*iPi, P 2 ( fzfc. L s O^Pjs P 2 £l, Pi + P 2 = 

1) fc-r&ts /JMusoaus^a. d&tizm. 

[0 178] 
[13 2] 

/(x) =P l / 1 (x)+P 2 / I (x) 



[0179] f ( x ) <&P*9M kftWiam, *tl*tl 

<!K&xmztih» 

[0 180] 
[»3 3] 

[0181] 
[»34] 



[0 182] Mi*Mt<nk%. Pi + P 2 : 

[0 183] 
[R3 5] 



1 SrfflV^T, 



o- 2 =-(a -^) s 



1+ 



1+ 



[0 184] Lfttfot, P x 1fi0frkl$.T&eflEth 
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[0 18 5] t>U 
[t3 6] 

[0 18 6] tth-h, 
[137] 

[0187] co^mmtz^tLhts:hU . cf2«SM 

j±, o<p 1 <i»ix:ra^#ft'f-^<, icoit*^, 

/hfli$#*x -/ -JayUzm-o X^ht^lz, **** 
[0188] fulBx ^gfcfr£(i>£-f C H 5 7 fc* 

[0189] A^!i«£ [ cr * * ] „ tti^W^Sr [ a * 
*] 

[138] 

/-> 

[ 0 1 9 0 ] #J£ffM«-C'0^%fgt^{i. 
[13 9] 




[14 0] 

[0191] ifc, W^tRjlz&lf&ftWta^cDimzk 
[0192] WjOMS&tJ; 0^5&^7 * ;^Sr# 

t.Wk=l ; sk2^sm2^TOme (1, 2, 
3 , 4 ) % 0 ; miBWhcOk £3tfftf-& & . 

[0193] ZtHi, 4~0<7>'hffi®<7)d%. t h 

X>y ^W«T^h««COW38^TO(C^f 
& COT . ®6T#iV vx -y S/ffi^j <I t &X% h i}K 

km. h^fr<nftWLn^.mx%m. i m<7)WMj) i mz. o , x 



[0194] aMBflcOlftSH'X.y i^SfiWi. ffi«J&»-£ 

[0195] 
[14 1] 

d" 
dx*dy"- k 

[0196] zcD7$-^]s-?te. m&x&mzte-ft&L 
^mmm^*^-* xh h . 

[0197] m.n&fotlh<7)\i, nifil tfcli2<WC 
n= l<0B*Hi, (3/5x), (5/3y)(7)2 

[0 198] 
[14 2] 

[14 3] 

[0199] f^7/W«^(i £ft&<7>»#ftt* 

[0200] 
[14 4] 

[14 5] 

*y(t-Jhf{*,/)-f(i,J-i) 

[020 1 ] Ate, 9&rC3k%ti&. 
[14 6] 

A = + (Ay f f 1 ,0 = tan - 1 (Ay / Ax) 

[020 2] 2&*#fcLTte. x#|6K y^ffl-HV? 

[0203] 
[14 7] 

A 2 *(/,/) = Ax(i + 1,./)- 

[14 8] 

W-Jhfit'J H) h/(;,y-l)-2/(z;;) 
[0204] -€-<0H*», Tia<D5gifi£ffl^i> - £ *>"T* 

& o 

[14 9 ] 

Ax(;, /)=/(/ + !, /)-/(/ -1,/) 
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[»50] 



[R5 1J 



Ay(/, /)=/(/,/+i)-/(/,/-i) 

I 

A 2 x(z, /)= Ax(i, /)- Ax(z - 1, /) 



[0205] ltftffcJHtAfcLT. 
[^52] 

|Ax| h|Ay| 



[K53] 



max 



[^54] 



(Ax) 2 +(Ay) 2 



[0206] 7r7^rvit 

[R55] 

J dx 1 dy 2 

[»56] 



[0207] Z.<7)1&* 1 g3&th t 
[»5 7] 

h /(/-!,/)+ /0',y+i)l" 
+f(t,j-i)+f(ij) J. 



= /(/ + 1, /)+ /•(« - 1, /)+/(;, / + 1)+ /(;, / - 1)- 4 ■ /ft /) 



v 



■ 5[(?(»,7)©ie«¥i§Hi:}- /(/,/)] 



[R6 3] 

5. 





f "* 


0 


r 






-r 


Ac = 


-2 


0 


2 




0 0 


0 




<■* 


0 






,1 2 





[0208] ztibawmmt. mm^it^ma-t 
^■^commmit^m^thtz^z. 3x3^11* 

[02 09] 
[R58] 

^ s C 

Z> i? F 
[g // /J 

[02 10] 
[159] 
1. 

A= {a+B + C-Q-II-I) 1 +(A + D + G-C-j<--l) 2 f 

ime 0] 

2. 

A = |^ :-S + C-G-/f-/| + |^l+.D+G-C-F-/| 

[R6 1] 

3. 

A = E — min^.S, C, D,E,F,G,H, /} : W. 
[S62] 

4. 

A=|S + h|^l-/| + |C-G| 

[02 11]Jj[TO5, 6i±. x, y«2J|S]«»I 



[R64] 
6. 





<-\ 


0 


i\ 




-1 -1 -r 




-1 


0 


1 




000 




"I 


0 






,1 1 1, 



[0212] ffna5fcj;tX607 jjufimzn ot, 

[0213] &tz, x , y cr>mjT$wffi.M&<omim.z 
[0214] 2&m^H;*M«-& ^r^s-ryc^M. 

[ISC6 5] 





'0 


1 0^ 




r l 1 




v« = 


1 


-4 1 


v,= 


1 -8 


1 




1° 


1 0, 




J 1 





^•fflv^^ti-g., ztuzMtxh. mmmzm^-ti 

[0215] mmjy'y^Tya. fom&ftttmLx 
iz£ y) s mm y^-7°% hcDi,z3mt& £ fc tf^tre 

[0216] Ztl*mm-& 7 A !V?l±. 
[R6 6] 
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'0 


-1 


°1 




-1 


-r 




-1 


5 


-1 


v = 


-1 9 


-i 




,° 


-1 


°, 




,-1 -1 


-h 



[17 0] 



[0217] ztit><ny 4 }V9<nmz^ {SMttzimm 

[0218] 4 /l^{±. iSttl£#H ( ju , v ) M 

[0219] m&y 4 jifimmiitipm^Mz 1 

[»6 7] 

fc-f*Ci:*C % EI58 (a) tf)<fcd%\ 7 4 /U:?l?9fSt<?) 
[0220] J*?—y—X7 M&y 4)V$<F> 

[*68] 



si' 



[02 21] n= 1^B#(07 4 )V9W38t.<?ffi&W5 8 
( b ) {C^-T. 

[ o 2 2 2 ] ifcte. Sf«7 -f /whtojh-*. 
mmt nmz, m&y ov? t w-v-x? ovtw 
[0223] m&y 4 )vm^ v&x&atih. 

[»6 9] 



[0224] A'^-V-X7 4 /I^f±. tfsST'fcSft 
[0232] «7c|8J|Hi:H:. ^-^H#S*^HH«f 
[»7 4] 

G(//, v) = B{n, v]f{;i, v)+ N({i, v) 



i i- 



[0 22 5] itt4><^7-f;^M»*05 9«O (a) , 
(b) IZyfrt. A'^-7-X7>f/^^t 
(b ) n=l«MTM. 

[0226] ztt^mm. w&yjfrfiz* ximw* 

[0227] <B«fl7C>miB7 4 )\s?tiZi. OfflDESfi 
[0 2 28] H«&^-^-r&/S««0«K{i, i6 0t 

[0229] Jn^WSry -f X<7)*§1=!\ f ( x , 

y) , /Utn (x, y) N MM ■ ffimKgtSt* B 

[»7 1] 

g(x,>-)=ii-/(x,j)+«(x,^) 

[0 23 0] BjMKJ&tftf&£W\ 
[1C7 2] 

g(x, y) = £ \f(a, p)>(x, a, y, J3 )dad/3 + n{x, y) 
[0231]C££ S b ( x , a , y, 0 ) ti. ^U£* s 

[t7 3] 

- «, y - p)dadp + n(x, y) 



[«7 5] 



F( J u,v)=M(/z,v).G(//,v) 



[0233] U «Jfto{^«l't'TJb-3fe*&<0^-^ 

(t) . 0 (t) fcu TS-isaji^Bifc-rsfcs 

[iSC7 6] 
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g(x,y) = 0f{x-a{i),y-fi(t)}# 



k&£>, 

[0 2 34] ^tlSr7-Ux^-rtL{f % 
[1*77] 

G{p, v)= Fiji, v)J^exp[- /2ff{^Jf(0+ 



[178] 

ff(//, i/)= |^exp[- /2*fria(r)+ ^(r)}]tf 

[0235] x^ft£-»KaTlflW~£«llK "T* 

« (t) = at s /S (t) = 005k$. 
[17 9] 

b(u, v) - — ^— sin Ma fit 
xafJ. 

[0236] B ( At , f ) Wfti)*-iXi3 *), 

Xi>ft&L%^i8r0rlz^ F (a*, Ojt5£ffif4, 
HR8 0] 

VM ' ' B{ji,v) 

t~ttiM 

[02 3 7] L*»U B (/x, »/) <7)ffi#04fct40fc 
ig^B#{4, ^OfTOOHdWIWfctfa k>RSii&fIl^ 
^ < „ »C J A Xtm& t X ^ & k * £ ttlSJS k * s . 

[18 1] 

k&D. HRfc. B (xt, »/ ) «N (At, u) fcjfr^ 
ts ( At , i/ ) #liD££«ft * kM^# < . ( m\ \ 
m%m.J£ttcoffi®X'te) 0K:}£< &&£>£*[ U J AX 

iJSsastfci'jfiSffiiv^tB (At, w^okjJt 

<*4J:,ffl*mKN (ju. ^)/B(a«, v)<Ofc*« 

•c«B (At, v) xmnz'rth-nz, ■tcot.t^m^m 

(B (At, v ) =1 kf&) ikkri-^T. MIBIS 
[0 2 38] -ftth-h. l7^;^«hU, 

[as 2 ] 



[0239] $ & J A Xffifo iHfetZ^ tmcowc 
BBR£fflV*£t><OkU »3§3rft5£fl7CORKi: IX 

ia. *m 2 mwm*m>nz-t& hmm^z i k 

[$C8 3] 



ZZX\ S nn (At, */) (4s y>fX»7-'Jxx^b 
;WT*0s S„ (At, i/)ttfif»7-UxWWI/ 

[0240] y>TX#iRV»*£\ S nn = 060 

[R84] 



f + v 2 > rt* 



[0241] CKOittt, y^x*wi>^{l ¥ 

i^m^MMXM^Wt^o tzMz 1/B (At, 

[0242] 4ft. ie«iai!i>-XT-A<o^%— tt^oaiE 
(54^ y^co* k PfiM H ) SIMfW&M 

[0243] t!l5gv-XxA<Dft«£ X "3 . 

(x, y) kf&k. f (x, y) 

j&*«8B£*l4^$k«I*K ^-tf)fctf>. 
[R8 5] 

g(x,y) = i(x,y)-f{x,y) 

[0244] ^«0^i%— ttSr1fljE-ri>^A6^!4, Mi.(4 
Sfe<0— ^J8K«[<OlHftf (x, y) =c (—56) £® 
g c (x, y) ##£>futfc-r£fc. A7JH# 
g (x, y ) (CMLT. 

[^8 6] 

[0245] MM^&liyE^^, ft«M*^k^^ 

m^fe. ^(x, y ) = ( r d , s j ) 
(x ' , y ' ) = (Atj, »/ 4 ) T&Sk-rttW. i*L^> 

[iS(8 7] 

jc' = ax+by \~c 
y - dx + ey + f 

fcrftAU #(«aSriR^tA(4\ S^OffliE**"srtgT^ 

[0246] tztzL, f^^jmmzMLxa. f^trn 
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^±^*v^ k & ^< s HfjS/JiSfifrfe 
[0247] <M£lk LT3^KK£ffl^fc£JBE8B> 

s. wait. fficoig&. jrjm^msmmi&z 20m 
[02483 mmmEmmt vxa, nmmmm^ 

<OLED^^#*><a!5:Sia^HftS6fi : LEDTW'^«S 

[02491 063(1 mMJyTfrb&hffiWBt 
tgllOl^t. H6 3fc*JWt, 10 2iiAny 

RSi-*-£fctf>tf>RtttKl 0 3 ^ffltiTV^ft. ^nyy? 

yno2 coffiMifftizii. jS.ffimmm-7 ov-?-\ 
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Abstract of J P 5244508 (A) 

PURPOSE: To provide a method for improving the 
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|f5l<lI141B«0^. 

c i**]a i 6 3 x n-r^-« «> t < » 

[USE l 7 ] NOEmR^. TIBs^ : 

F (x) = -m* ( - x /m)P x <0 CO*^ 

F (x) = m* (x/m)p x >=0 

( ^cfjgp ft 0 < p < 1 <y)iGlllF*I-C-»K§tL. m^f Jl 
fStt^iEH : -m<=x<=m£^rr) k LT^«§tL 
■Sff^Hl 1 1 KEttO^fi*. 

[ff^if 1 8 ] HBEHfttf. TIEjft : 
F (x) = -m* ( - x/m)^ 2 -m < = x< -c 

corns 

F (x) = -m* (c/m)^ 2 * ( - x/c)^ 1 -c < = x<0 CO 

F (x) = m* (c/m)p2 * (x/c)^ 0< = x<c £0^ 
PI 

F (x) = m* (x/m)^ 2 c< = x< = mCO 

^ P 2(iO<p2<lcy)iEH^T'jlfKS^. X 

C0K#jS : 0 < c <m SrS LT*3 i.Xfm\iWMUcom 
SI: -m< = x<=m^-r) k kT^«Sa^lf*Jl 
l 6ta«^^rffi. 

1 8(ciB»<0*t6. 

[If 2 0 ] flulB-7 y \z°>7Mm\ SWE*IP*t= 

#«wrk u ^-ch t tttjEii § fit* t= s <t ") mmm 
msM i * 1 1 ciBSJfXfris. 

[IS*^2 1 3 »I$^H«cD*M -y ?l/>y 

^\ mmztitimm* . mfecommmkiz tt^^x. 
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BtMOtfJ^'-f ts-y?i/y i?t,Z-? >y t° y 6 £ k £ 
X o T«6 SflS 1 ~ 2 0 cr>\, i-ftifr — oiztdM 

X 0 i& v ^JWMEU ^Wzfc ft 4 S^BfR h ¥ 
K, JftW3IBMttfi***rf £> iz~?tix&« lzm&tZ> xu 

-vk^-$h'}?ti:< t i>-m$mBTmmi,zM±'t& 

>y try ^jRt i T , «riB^H*tf>iB*«£ 
r;i/3" u Xi,^ affl-f s £ k x -?T*aa$*utH«& 

msm2 3 ] m^mm&u^Mziditfmmmi^ 

<HlB«Srlt3t'^S#S, m&cD^MztslfhjZimm 
k &tf)fflv \\s-xjuzt5 ft h 3fi<KB#£— § * 6 ffif fc . 

{za (t h mmmi&^-^x o «v «««y ^/Hcfc (t 
h mmmztm l r mmmm:-^, $ s mat #s 

^WMm&xm&LxR&r . if**! 2 

2tmco!£W.<, 

{wmm 2 4 ] mjmmim-t&mttmK vm 

{ it *if 2 5 ] m y r u y 2 t«£ t 
^y^-yrv yyimz^-tz £ k Sr^ti&k u itria 

i£«7 -f ;W^*szi<X7t^X^S:^rtl.|»*3l2 4 IB 

[fi*ii 2 6 ] ^m%ms.^Mza ft & **bjmb 
mtm^m^m^t ^■t^mm^mm.^ 1 * -y b & 

h&^-m, m*mvk<7>* -vv^o *>^> — 0 
»w I* -y h * . friBH^^^m-rs - k (C «t o-CfW 

•T-l>¥STS)-oT. fTIB^{i> mrlB^m^ife^SfflS 



xu-rz^-tz, >y+£< t h-mimmxnmizm 

^tLtzmmimtn^Mmz^oms., mi$. 

[ft*iB2 7 3 tmmstftx 9 y -y £B&#tt»st 
k (c x ~>xmx9 v->wzimznx 

* § 6 (zM fl^s 2 2 4 fcli 2 6 IBifcOiSB. 
[000 1 ] *5MHIi-r -f : J?Jl>mWcoa yhyXhi 

mttimm^yhyxb^mmtrnm^zm-t^. $6 

tcpL<{i3yb° A -^ ftfftftflt 1$ yxT-^SfcJiny 

&Bisx^j±i l zm^&m?.cojjmi l zm-t&. 
[0002] mm&jjMmt vcommm^mmt^i-f 

ixS * U ^'-t/WB^^JiSv ^ISOieHrt-r-a^ cOfi$ « 
fis ^^^^s-y^ixyi/tc^-p-caiBfJrb^t 

x^yy-tf&imzn-ox^&tz {xy\?x.~srmt 
yb° ^-^mmmmm^m^-h^ a 
■c*ft ) tz*ai -g. . z\ cox 0 %mwxt 6 1 ^titcfm 

fl^B#. yt°jL—?%MWWt&Mi$. 
fiK 7 4 n — K 3 f — (7^f b -y ^ X±T-^IS 

iis ) ±*>*fctd*&«^T«>7 f 4 ^7V-f u-y 
±{c^Sii »|giJ^±«BJ*H«3 y h 7^ b f4S 
fc^-^ t- 5 -y ^ y y v'{c>!t LT b y- H ^n^rfttLff^ 

X\^Z>i9& (Mm? 4 ;PA<?)Ji^-500 : 1 ) XK)'hZ 
< , a^O^M*#T<7) CRTX^ >J— y<0*& 100 : 1 

i 0^h§^) M^-H^-7{iSSBtc^yxlf OT<0k*J 
0T-a5«». 

i ) mm±Mkcoh h m^^vco^y^ t- s -y ^ y y 
y>'f dffl t d « xkw&wa ts-v^y/ ycv -y r 

$ft&*ily ^3Vl5Xl*Wft<. ^<^ff^ 
BHWtcoi \tmi t # v^y<;u i 0 /h^ v ^TH^# 

ii ) —73^ u y't^y^f -t s -y ^ y yy'o— i5fifw 

v^Vy^z^-zT^fih^, t\ 

^mvxy<n±m^v^M±^x n t «£vo aitjy 

[0003] ^c0^-&. miR$tlfz¥4 tS-y^ y^./k 
rttc y^WSr^f-o X v> S ^-fi fec7)H^B*!^*ft^§flT' 
§&rjyb5Xh£ J ?-;t6:f-u — fficoB««±«ft^ 
%g—%:cr>X'-£< ziyY^X h~&LXyjk£fi& . 



(5) 



#§PP5-2 44 5 0 8 



[0 0 04] nyb^-^MWiaKvX'f At L<(i 

[0005] L#>L«7)J: ?&»0#itiS*rCf4. 
SjI^/K/I/ BAN"* BSH* 1 * v * . L£ 
T s fflmiiWbX SBBK*- <0H« £ 7 * /PA & L < {4 

[0006] i^mix-f : J^^mMn^^mxit 

&<%atztlX^& (Maack I., Neitzel U. , "Optimize 
d Image Processing for Routine Digital Radiograph 
y" .Proceedings International Symposium CAR ' 91, 
109 H Springer Verlag #JS) . 

[0007] nakono^siR ■ «ft* it^ryr^ 

flttD 5 f-f - jl — K #i4 & V •> HZ #>#>;b fe-f . :£ 

[0008] zcom^mss.-th^nbiz, -mmz%\t> 

itX^h^^TViy^—y^^yf (unsharp mask 
ing ) s ®j£6DbXb^7AW-ft: (adaptive histogr 
am equalisation ) &£X/Ztlh<D--mmijm<7)%<<7) 
3cBcr> Xd^< comfrtfftZtlX £ fzi)K -f<X » T 
-f7T^b ( arti fact ) t m$tl& Zf—A b H ^B**, 

fohmmfo *) . - tumnmmizte if h 'W/ms.mco 
n$tzmm&m*i&Fr?& ttzimco 

^*ffffl{4S3Kfc**KI&«S^TV^ (RehmK., D 
alias W. I., "Artifact Suppression in Digital Ches 
t Radiographs Enhanced With Adaptive Histogram Equ 
alization " , Proceedings of SPIE, 1092 ^ Medical 

Imaging III , 294-296W 1989^ International So 
ciety for Optical Engineering, Bellingham; Oestman 
n I.W., Prokop M. , Schaefer CM. , Galanski M. , "A 
rtifacts in Digital Storage Phospor Radiography 

" , Proceedings International Symposium CAR '91 1 
25 ~M Springer Verlag; Bick U. , Wiesmann W. , Lenz 
en H., FiebichM. , Von Lengerke H.-J., Peters P.E., 

"Utilizing digital luminescence radiography in pe 
diatric radiology: a report of initial experience 
s" , Electromedica, 59 m 1# 30K 199130 . 

[0009] ztit><7>3 yh?x hm-itmcom^mm 



14. JBBiaftn— ^si-^u— ^SttdiktSv^Tv^i 

tiK ^\s-?co±%^cr>mmzwmthi><7)T$>h . 

imttstiz . mmttMmmwiWMmmiz^v- 
fzm<mxi>) T»srofc/urrjBi6. 

yynzM LTllfSt" IBM'*?* —9 lz 

ov vc < cotmm maz i^tr *u4a**Rnrc* 

S£fcjWIS3*VCV*4 (Prokop M., Schaefer C. , Oes 
tmann J.W. , Meschede A., Reichelt S., Galanski M. , 
"Optimal Parameters for unsharp Mask Filtering i 
n Digital Chest Radiographs" . Proceedings Intern 
ational Symposium CAR ' 91 149-154K Springer Ve 
rlag; Prokop M. , Galanski M. , Oestmann J.W., Von F 
alkenhausen U. , Rosenthal H. , Reimer P., Nischelsk 
y J. , Reichelt S. "Storage Phosphor rersusScreen-F 
ilm Radiography: Effect of varying Exposure Parame 
ters and Unsharp Mask Filtering on the Detectabili 
ty of Cortical Bone Detects " , Radiology, 1773S 
1 ^ 109-113 H 1990^10^ ) . 

[ooio] t< vfivmmmnftwx*. ^mmmm 

ft# (multiresolution computation) COfftlfrtffc&^i&S 
tihb^hfoh* ZCDWkl&lZktUi. &%LCD*>yh<7M 

m^7*-7tfmm$tiizJKm.mcDTJ s>*/hm xt= 

<X<0S3iW^r*>^. Burt P.J., "Fa 
st Filter Trans forms for Image Processing" . Com 
puter Graphics and Image Processing, 16 % 20-51 
M 198W; Crowley J.L., Stern R.M. , "Fast Compu 
tation of theDifference of Low-Pass Transform" . I 
EEE Trans, on Pattern Analysis andMachine Intellig 
ence, 6^ 2# 1984^3 M lz95M£tlX . fl&co 
fliMWlSfcltii, Mallat S.G., "A Theory f 
or Multiresolution Signal Decomposition: The Wavel 
et Representation " , IEEE Trans, on Pattern Analy 
sis and Machine Intelligence, 11^ 7 # 1989^ 7 
J^; Ebrahimi T., KuntM. , "Image compression by Ga 
bor Expansion" , Optical Engineering, 30 ^ 1 ^ 
873-880 M mWlBifihh, 

[ooin cwmnmimimffico^mg 

WH&Wi&COEfflX'$>r>iz (ArbeiterJ. H. , "Multidime 
nsional video image processing architecture" , Opt 
ical Engineering, 25 & N 7 ^t. 875-880 M 1986^F 
7R; AdelsonE.H., Simoncell; E. tS&lf Hingoran; 
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R. , "Orthogonal pyramid transuforms for image codi 
ng " , Proceedings of SPIE, 845® 50-58 M 1987 
International society for Optical Engineering, 
Bellingham ) „ 

cooi2] mammzte, ^mmmm®.*.?** yf- 
-~> a y , rn&mm s *j J; tm#iwjs»»jesf i= i & 7 

4 )V? (Lifshit, L.M., PizerS.M. , "A Mu 
ltiresolution Hierarchical Approach to Image Segme 
ntation Based on IntensityExtrema" , IEEE Trans on 

Pattern Analysis and Machine Intelligence 12® 6 

-f\ 529-540 M 1990^- 6R; Szeliski R. , "Fast Su 
rface InterpolationUsing Hierarchical Basis Functi 
ons " , IEEE Trans, on Pattern Analysisand Machine 

Intelligence, 12® 6 # 513-528 H 1990^6^; 

RanganathS., "Image Filtering Using Multiresoluti 
on Representations", IEEE Transon Pattern Analysis 

and Machine Intelligence, 13® 5 ^ 426-440 M 

199W 5M ) . 

[0013] *%M7)£ftBtm^ yj+^vtvyv 

m^V^V^Vy^zhtz^X 

&C\bXhh» 

[0014] ajc^oflfioBWJi. i&mmcommcDay 
h ^x h £ i b % < + S y 9 V y : J*'h% < 

tfX% tztz i k T*> & . 

[0015] *^Bjc7)ffic7)aa<j{i, »jm>mh*>x*;S § 
-r^xdf y^t fcttSSiew t x h /7Aff»ra 

? b £ &J&t * i i: sfir < 3 y r- 5 X h £ 4 £ S 

[0016] *%iW>% ^tffieoaWtts SfMlfrfc* 
t < k&^yb^Xr- £5MrT£;«SiH5« 

tSiktibi. *%BJ<50$^tCfB!<7)eW(i^$ (Shar 
pness ) £i|!*£*&#63yh5Xb£»<^.&#& 

[0017] 4>teJMcOBK»4JJEW»t*t 

[0018] *^<0»BB#fefi. H*{fi?»K"St§ 

jsmsritra-r s £ 1 1 «fc y Hi^r? yb7Xh&a 
-fw s^rffittjv . ( a ) imb* y ss-j-Mpme . 

^mv-^JVco^bX K) i><3\< V ^m&bz L- ^Mz& 
^WB. b izttM ( b ) |frlB«0J«BB«^H«« £ 51 

i>^< fc is-r>coimwx*wmi l zMtt&mm.wm 



£#t , ( c ) mi bs#h« m bodies 
titzmmmzmm$LT)v ^ y xa. zmmt hc\t\,zz. 
^xmziowmztifcM&zmn-f&ztfr z 
emmi&TiV'd u xa nhbz turns.* t xM^mu 
b mmmzmmztitz* he* u ^;m*4 ^ 

<b-r * i k *= «t -5 t . ±ieg wjijftt s 

[ooi9] ttrie^wji. lufB* u ; ji~ii>mi&e>%mm 
m&m?.&mm&cottiB~t & mmmz x ~> Ttgjnstt 
tzm^mm^<m^hm%m.(ry^m%mzm t < 

wmm:msfjMMz x^xxv i>i-)vmmz-$tz 

Zcom%=. i ) *ib»>tWWW«^:H«*W& 
ffl^-fe'nic^rD ; ii) fe^^SBJfflH^SSil^C^ 

-otjkAk-r^T (#990%) **:rrssisw8i!sgRK 
xjyttxD&timmb Lxgmzti. mmmmm^m 

>J ^/l^fflSWiHR^J^KtclWIfSti. ; iii ) M^O 

Wttz2>& b ^^titz-t^xcomm^mmzm^rf 
h m^mmmmx-t) ^~^tih^x\ 
mmmmm^h^mmmm^m^znmL ■, iv) iuibbj 

m^bwm-thmmmmiz^^itttii - 1 a* 

v ) #^B«rt^H*<^R**iJ !,i 5:< b i>ffi*)ML*M 

ifhtzMz^ 1-4 ^x b-^-yy u y^ap^i^TW* 
sfLftB*<?Dirr* o ; vi ) liriESiajaasdRWW^^ 

[0020] ®ffl§ixfcMm«. Wm&OftR&hM 

mmxm^ti^^miEmmm^mm-^m b txtm 

mcoms* y i?i-)vmm.<mm<F>mzmL< 
m®%mmz£-?x-&mt$> . 

[002 1 ] *5MB<oaW«. § kfcJTFOXf-'/T^ 
ZZbWCZh. -tftifh* a) frlE*y^'^Hli« 

Sr. ztimmt&zbizx^x, ^mmm.'&v^Mz 

nnMrnsmzftML. mm^mizx^x. mtmm 
Mfo*$?t i -fe -v bcom^mwtco a — ?co* y e^-^ 
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9 s flufB^lte. lw!B^}&«ifc«iffl§ftfti§i^ *U 

£-tft £ fc U b) ^*f5l^ffl^ift^TS tco 
ftT#*tcS^^ftxo-r£^ft^3:< fc — ?0 

T«ffi9Bffl3*ifcM**tfltt"* , x-f- -yT^&ftM 

[0022] ^<onfHHS(z imr. mfB^ossit 
me* u ^;m«i*i^)#n*<o«[Ji. mm-&wm® 

6tut«*HJWB*co^*-b -y N**s PI t 

4HW $ ft ft *£■ s lift B^raiKH^fiM^Sfflffl*^ 
^tfsyyN 0 ? hB#Ri:5£ifi$tu »«IB 
h K^M yco*t»0:£BS^nffi£3rU ^ 

tomb ^titzm*mmmmmwL<7)Z. a tc4HW$*ifc*s 

5&*-fe'n(c^r9 Mi) &&»9>4 iiufB-fr-te >y h^mflB^U 
¥1-)Vm<$tC03L K .yM y A- U -f t> , MIB^j" 
Ui^/l'K.XW yftco&^ftBSRfc-ovvr, £>fe»9>ft 
BulBl^-fe >y b F*J£ . SHISB^iiuIBB* t MM? ft & 

•fe h kuw § mt b**bhw . pi t^isua 
a^-t^t ($«j9o%) z&tsgmmwmvxj y*<oa 

yn? hffl&b LXimtZti. iv) BSP B lJtMfflt& 

^x-K—Ttcriji/T y/mm^mz $> h ^smajs^ 

■T^T fcBBW SMfBSr^^^^"C^/N*-$iiS J; 
at, * & »9> ft mifty* >y h j^ 5 * ft SlfflJBaa&WJ* 
fc*tJEU v) lu=B^7"-fe-y ho— ofcRBS-$-4#ffiia 

CO friBSPBlSIMSSc^iS^mrlB^SlOiiS £ & * fc 
ft ft. 

[0023] *^0j{±§^t. mfmmm.^z^x , 

fcHR*?-* . *o^3ftfc IWiRii; v- H ^ t°-IB#i. 
i/Hit-^-cofari/'fct^t bsj-cs ft £ a 



fc&ft. 

[0024] «-?B«lS*<ft If <£<i!K4*:i±«?-e< 
?y 3 yT1i BuiBWfSlftii. flilIK;W:B(t 

titzmmmm^T)v d u xMz£~oxa-M.zti&„ ok 

[0025] ±IBO«-e < t^yg yiiHtf»0«r?fl^ 

timmibtL&mziiXtmffi&E&s mmmmmm. m 
mmwmy < )Vj±±tz\mmm&wtx7 v-ycox. 
v\ 

[0026] liijl-t? is a yiii/-fWgl4fciiI 
IftefMS^ i: O J; d \- Y 3 f-IBIi^B-C'* o T 

[0027] *5KBoii««yis-ffi{4. mm^m^tx 
^v-ycommmtnL^w^i-<im^-t^coi>z. ~> 

VOVcr&^mzhiz-oX , -r i??J\sM1&cr>ay 
Xh£?m't&<?>£BWt LXm%%titz<7)Xfo%>. % 

[0028] *i&B^c7)^g^M-rft^o?hop«a*ii 
^wm^yhyxbmtmcom^&mmimfizm 

[0029] *m i M<7)mw.<Dffi%yv <y ^asria 1 

f . B««-?-<^Kli4> CCD -fey-9". tr^yf^i, 

7°9 y^'f ft i i: fc<k -5Tf v^;Him2rfI^ < U 

U 5^-;UWR 2 t Df (fflHft 1 8 — 1 2 CO b' >y h *«|ffl 
«ft * ft BM<7»J t , A / D Sr ffl V^T * 
L . BfrMfc «t 9 * coB*tt i * =E U fc— BSttKEfli 

muzx-oxmmmzmtzti. mzmt^titcmi^4 
tfayh^xh mmzz-oxmmmzfeTE-t&mTF^ v 

fy/t? y 3 y 5 icSMSti. -ecois*. «sa$it3t 
am 6 ^ h V < itWciZJ \- K a f 

9 4 /tii^ a7W^'J-y( CRT) ±lz 

wrtg5rB« ^#ftB«m*fi7 jciwr* ft <t a t$ix 

ft t # fc, ^T4B«flB8* I «a<0^rtS"TffiftSn 

ft. 

[0030] ifttffi* < i co»* Lv^SHB«*s 
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KSHftWRa*. «i'ffiSb£:? c l/-Ji£ . XWM1 0*> 

Lxm$m& * wt tfz x.mzmmt& ztt,z£~o 

fcioTffiitiSfi*. <I<7)f^M{4. IflSSSflMHSl 
SKi-j-cijaarifiifciidiSiiS. -<XfeBB»(=J: 

t J9)V$8ffli 1 9 (ci->T 1 2 t*>y t-O-r 4 S^/Hi* 
-JvMB* k D?«'tL Sf^'^ /WBft 2 # . Sfrfc-fe ? is a 
[003 1 ] WMMki/^'f^ 3 l4H-o<J9£ggfrJK1S 

wc. ^'jyt;«t2a, -n&cvWMmmz-ftMz 
x<r)$>wm®& u ovx-% v : Ji-)ummmzi?j£-t & w 
[0032] tn«xf 7r«a, n^mm3 1 - 

xiwm~? -v t-y^mmzi. ^xmiEztih . mmm 

Jg-te ? >- 3 y 3 4 H£ V vf „ 3fi-ft,Stt^ll#4 ^-tf-fiL-r 

k k t ^j5?®«L-^T'Ha§^^ . 
[00 33] ^«BJM^5*F4 U^HBfi»«£EB4 fcjj* 
•f. ^ U i/^tVMWfifeM? 4 )V? 4 1 \,Z X -o X 7 4 )V 

*s*u Mz 2 m-y-yy-yy uy^ztii. c\mmh 
»9> & ffi^HS^fia^*3 v>t co^-im-r h <r k fc «t -5 "C 
[0034] *imiffl^K;Hc:t3tt-I.BJ«ffllB«bo{4, 

sk ^j^iicS: 2 fg l , mmm^zmmmm^yf u ^;nf 
« 2 < <r k t i o t . ismmm wovmwmm 

[0035] ^owiwreB^ -oa 
tjzfc^&i&comzwAL. mz^cowzk^tuzm®. 

ZiSM? 4 ;lsfX°m&Z J £& (Convolve) . C\CDB\£Mte 

MMM4 3Xlrffrn&. 

[0036JH tWifttf. * i> *y*t)i-mm 2 <oftb 0 ^ 

it I.. 



[0037] —mcowfflmmi < i = o-l- 1 ) k^ 

fi . fttt(OBJ«bo <0 1 / 2 <J0a<Off k n t *» -5 T V -»=5r 

Si)SR«M^ta^lHiBffi«<0^*iO 1 /2 te3i*Ts 
i^ijkff^i4^-^ h»(:J:o-r 1 /2T* 

s. mmwM<r>^ -mm&zi3 1 - - 

^x{4^ u ^-^/Hifli 2 (n^^m^^t — mwwJSfx 
[ o o 3 8 ] Mt L^mmm<r>&®7 4 iv^coy 4 iv 

m^-f^xcox^r-^xi&my ^M4K 4r . ••• 

4 1' ' ' fcttLTfflV^fl*. 4SrSHmt7t-r^T 
2' •■•4 2' ' ' Mzm^btlZ. 4<7)ffim&±MW 

[0039] ttm^&mmmmzmsiz^-t. wm 

mm l fcaiBHS 5 1 i o x Z cn%<r>±% § <o 2 f§£ 

^§^B«4^#feix^. »13i5 1 s 5 1 ' -5 
1' ' ' \$ttffi*i9>>syX'm~?tz<ntnii><r)Xfo 

[0040] ffimftmmtm$ti&±MicoftMmtWM 

BStSO§^^f¥iffl ; ^P (: I^i4. "Multiresolution Image 
Processing and Analysis " , A. RosenfeldHS. Spri 
ngerSeries in Information Sciences, Springer Verla 

g 1 0 ~ 1 4 mztm^tix v \& . a^wsataex ^ u 

-yoMfM«S!^aB L«B£«iWfe I F H «f42 
OS : -¥-Srft*-&a i *^«4PfiJ^.i.iktJ;^TK»-r^ 

i k 4 . § ^izmm^z^Ki, t . wimiffl<&& 

OHf v ^I«{4 9 BBS v -fflturrtiif— mzm%1Zfrt>X%> 

h. z.co3m*mffcf&fcM^ mzMi^htifzmm 
mi. Ttz<vmtmmcodib<D — ^zx-oxmsE^ti 
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[0041] iw^acon i commmzmnz^k-t. 

mbBMtl&Mkto < * 4 1 < 14 2 *> t < (4 3 O ) ®pg 
«MJE-TI>H«t»m^5 3" -5 3' ' ' £2tfHfcJ: 

3" •••52' ' ' <,z i. ^xm.M^tifz~? -y \zvy\z 4 
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(57) Abstract 

A method and apparatus for the fast detection of spiculated lesions 
(1608) in a digital mammogram, the method for use in a computer aided 
diagnosis system for assisting a radiologist in identifying and recognizing 
the speculations among a multiplicity of lines corresponding to standard 
fibrous breast tissue. A line and direction image is created from a digital 
mammogram, and a region of potential intersection for substantially 
every pixel in the digital mammogram image is determined. The region 
of potential intersection for each pixel is a predetermined pattern, such 
as a high aspect ratio rectangle or trapezoid, positioned around the pixel 
and rotated in a direction corresponding to direction information for 
that pixel. The regions of potential intersection are accumulated among 
the pixels to produce a cumulative array (1606), and information in the 
cumulative array is processed for identifying spiculations in the digital 
mammogram. 
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METHOD AND APPARATUS FOR FAST DETECTION OF SPICULATED 
LESIONS IN DIGITAL MAMMOGRAMS 



FIELD OF THE INVENTION 

The present invention relates to the field of 
computer aided diagnosis of abnormal lesions in medical 
10 images. In particular, the invention relates to a fast 

algorithm for detecting spiculated or stellar lesions in a 
digital mammogram to assist in the detection of malignant 
breast cancer tumors at an early stage in their development . 

15 BACKGROUND OF THE INVENTION 

Breast cancer in women is a serious health problem, 
the American Cancer Society currently estimating that over 
180,000 U.S. women are diagnosed with breast cancer each 
year. Breast cancer is the second major cause of cancer 

20 death among women, the American Cancer Society also 

estimating that breast cancer causes the death of over 44,000 
U.S. women each year. While at present there is no means for 
preventing breast cancer, early detection of the disease 
prolongs life expectancy and decreases the likelihood of the 

25 need for a total mastectomy. Mammography using x-rays is 
currently the most common method of detecting and analyzing 
breast lesions. 

The detection of spiculated, or stellar-shaped, 
lesions ( "spiculations" ) in mammograms is of particular 

3 0 importance because a spiculated breast tumor has a relatively 
high probability of being malignant. While it is important 
to detect the spiculated lesions as early as possible, i.e. 
when they are as small as possible, practical considerations 
can make this difficult. In particular, a typical mammogram 

35 may contain myriads of lines corresponding to fibrous breast 
tissue, and the trained, focused eye of a radiologist is 
needed to detect small spiculated lesions among these lines. 
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Moreover, a typical radiologist may be required to examine 
hundreds of mammograms per day, leading to the possibility of 
a missed diagnosis due to human error. 

Accordingly, the need has arisen for a computer- 
5 assisted diagnosis (CAD) system for assisting in the 

detection of abnormal lesions, including spiculations , in 
medical images. The desired CAD system digitizes x-ray 
mammograms to produce a digital mammogram, and performs 
numerical image processing algorithms on the digital 

10 mammogram. The output of the CAD system is a highlighted 
display which directs the attention of the radiologist to 
suspicious portions of the x-ray mammogram. The desired 
characteristics of a spiculation-detecting CAD system are 
high speed (requiring less processing time) , high precision 

15 (the ability to detect subtle spiculations) , and high 

accuracy (the ability to avoid false positives and missed 
spiculations) . It may also be desired that the spiculation- 
detecting CAD system also be used as a mass-detecting and 
mass -classifying CAD system, and that the CAD system be 

20 capable of using spiculation information in conjunction with 
mass information for identifying suspicious masses in the 
digital mammogram and directing the attention of the 
radiologist to both the spiculations and the suspicious 
masses . 

25 One method for detecting spiculations in digital 

mammograms, proposed by Kegelmeyer et al and referred to as 
the "Alignment of Local Oriented Edges" (ALOE) algorithm, is 
described in Kegelmeyer, "Computer-aided Mammographic 
Screening for Spiculated Lesions," Radiology 191:331-337 

30 (1994) . The ALOE method first calculates local gradients in 
a digitized mammogram. For each "candidate point" in the 
image, a predetermined window around that point is selected, 
the window size being some fraction of the overall image 
size. An "ALOE signal" for each candidate point then is 

35 calculated based on information in the surrounding window, 
the ALOE signal being defined as the standard deviation of a 
histogram of the gradient directions of all pixels in the 
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window. The next candidate point, offset from the previous 
candidate point by a distance corresponding to the desired 
resolution of the search, is then considered. 

Keeping in mind that a spiculation is a roughly 
5 symmetric set of lines radiating from a central point or 
region, a histogram of gradient directions will tend to be a 
flat distribution from 0 to 360 degrees if a spiculated 
region is centered around the candidate point. Thus, because 
the ALOE signal is the standard deviation of the histogram, 

10 the ALOE signal will be lower for those candidate points 

which are at the centers of spiculations, and will be higher 
for those candidate points which are not at the centers of 
spiculations. After the ALOE signal is calculated for all 
candidate points in the image, local minima in a plot of the 

15 ALOE signals are used as a basis for identifying 
spiculations . 

The ALOE algorithm has several disadvantages. The 
primary disadvantage is that, in addition to spiculations, 
many unwanted background objects can also produce a small 

20 ALOE signal. For example, a point that is surrounded by a 
circle, such as the border of a circumscribed mass, also 
produces gradients in all directions, and therefore will 
produce a local minimum ALOE signal. A false positive may 
result. Furthermore, a typical spiculation in an actual 

25 mammogram will not have lines radiating in every direction, 
but rather will have lines radiating in several discrete 
directions in rough symmetry about the center. Thus, because 
every direction may not be present in the histogram of 
gradient angles around the center of a spiculation, the 

30 standard deviation of the histogram may still be quite large, 
resulting in a larger ALOE signal. This spiculation may then 
be missed. Thus, the ALOE algorithm represents serious 
practical problems because it may yield false positives and 
may also miss certain spiculations. 

35 The ALOE algorithm is representative of a class of 

"backward direction" spiculation detection algorithms. By 
"backward direction" it is meant that a "candidate point" is 
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incrementally moved across the image by a distance 
corresponding to the desired resolution of the spiculation 
search. At each candidate point, a set of "window 
computations" for a window of pixels surrounding the 
5 candidate point is performed, and a metric corresponding to 
the presence and/or strength of a spiculation centered on the 
candidate point is computed. Thus, for example, the ALOE 
algorithm computes the "ALOE signal" for each candidate 
point, and then moves on to the next candidate point. 

10 As a general observation, "backward direction" 

algorithms are computationally intensive. This is because, 
for an image size of N x N, there will generally need to be 
on the order of K (bN) 2 computations, where K is the number of 
window computations for each candidate point, and where b is 

15 the reciprocal of the number of image pixels between each 
candidate point. Because the number K is often proportional 
to the square or cube of the window size, the computational 
intensity of "backward direction" approaches can easily get 
out of hand . 

20 A second method for detecting spiculations in 

digital mammograms, proposed by Karssemeijer et al . , is 
described in Karssemeijer, "Recognition of Stellate Lesions 
in Digital Mammograms," Digital Mammography: Proceedings of 
the 2nd International Workshop on Digital Mammography , York, 

25 England, 10-12 July 1994 (Elsevier Science 1994) , and 

"Detection of Stellate Distortions in Mammograms using Scale 
Space Operators , " Information Processing in Medical Imaging 
(Bizais et al . , eds . , Kluwer Academic Publishers 1995). Like 
the ALOE algorithm, the Karssemeijer approach is also a 

30 "backward direction" spiculation detection algorithm. 

In the Karssemeijer algorithm, a "line image" and a 
"direction image" is first formed from the digital mammogram. 
As is known in the art, a line image contains line 
information for each pixel in the digital mammogram, while a 

35 direction image contains direction information for each pixel 
in the line image. The most basic form of line image, used 
in the Karssemeijer algorithm, contains line information 
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which is a "1" if the pixel is located along a line and a "0" 
otherwise. The most basic form of direction image, also used 
in the Karssemeijer algorithm, contains direction information 
which, for those pixels having a "1" in the line image, 
5 equals the approximate angle of a tangent to the line passing 
through the pixel . 

Consistent with its "backward direction" character, 
the Karssemeijer algorithm then considers a window of pixels, 
chosen to be an annulus, around a candidate point in the line 

10 image, and then computes a metric associated with that 
candidate point. This procedure is repeated for each 
candidate point. The metric for each candidate point is 
calculated by counting the number of pixels in the annular 
window which are contained along lines which point 

15 approximately to the center of the window. Local maxima in a 
plot of the metrics are used to identify spiculations in the 
image . 

The Karssemeijer algorithm has several 
disadvantages. The primary disadvantage is computational 

20 intensity due to the "backward direction" character of the 

algorithm. Because each line image pixel actually appears in 
many sequential windows corresponding to successive candidate 
points, the calculations are repeated for each pixel many 
times, and the algorithm is very time consuming. As 

25 discussed previously, speed is a key factor in a CAD system 
for assisting a radiologist. A CAD device which slows down 
the radiologist, who must often view and analyze hundreds of 
mammograms per day, is undesirable. 

Further, the approach used by Karssemeijer in the 

30 detection of lines for generating the line image is based on 
using Gabor filters in the frequency domain and performing a 
Fast Fourier Transform (FFT) on the mammogram image. The 
image pixels are then multiplied by the transformed Gabor 
kernel elements, and the inverse FFT is then used to obtain 

35 the enhanced line image in the spatial domain. This approach 
requires input images which have dimensions that are a power 
of two, which is required by the FFT. Thus, after digitizing 

- 5 - 



WO 99/63478 



PCT/US98/11729 



the mammogram, the Karssemeijer approach requires the 
digitized image to be "padded out" to the nearest power of 
two. The cost of this approach is higher memory requirements 
in the computer and a larger computation time. 
5 Accordingly, it is an object of the present 

invention to provide a fast computer-assisted diagnosis (CAD) 
system for assisting in the identification of spiculated 
lesions in digital mammograms, the CAD system being capable 
of producing an output which directs attention to spiculated 

10 lesions in the x-ray mammogram for increasing the speed and 
accuracy of x-ray mammogram analysis. 

It is a further object of the present invention to 
provide a fast CAD system for detecting spiculated lesions 
which produces fewer false positives and fewer missed 

15 spiculations, while also being capable of detecting smaller 
spiculations . 

It is still a further object of the present 
invention to provide a fast CAD system which detects 
spiculated lesions in a manner fast enough to permit use of 

20 the CAD system in a clinical radiology environment. 

It is still a further object of the present 
invention to provide a fast spiculation-detecting CAD system 
which is also capable of detecting and classifying masses in 
a digital mammogram, the CAD system being capable of using 

25 spiculation information in conjunction with mass information 
for identifying suspicious masses in the digital mammogram. 



SUMMARY OF THE INVENTION 

These and other objects of the present invention 
30 are provided for by an improved CAD system for detecting 
spiculated lesions in a digital mammogram image using a 
"forward direction" detection algorithm, as opposed to a 
"backward direction" detection algorithm, for improving the 
speed, accuracy, and precision of results. A CAD system 
35 according to the present invention employs a fast method for 
detecting spiculations in the digital mammogram image, the 
method including the steps of determining a region of 
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potential intersection for a plurality of pixels using line 
information and direction information related to the pixel, 
accumulating the regions of potential intersection to produce 
a cumulative array, and using information derived from the 
5 cumulative array, such as the positions and strengths of 
local maxima in the cumulative array, for identifying the 
spiculations in the digital mammogram image. In one 
embodiment of the invention, the region of potential 
intersection for every pixel in the digital mammogram image 

10 is determined and accumulated into the cumulative array. 

The line information and direction information are 
obtained by generating a line image and a direction image 
corresponding to the digital mammogram image. The region of 
potential intersection corresponding to a pixel is found by 

15 determining, according to the line information related to the 
pixel, whether the pixel is located along a line, and if the 
pixel is located along a line, selecting a region centered on 
the pixel corresponding to a predetermined pattern, the 
predetermined pattern being rotated by an amount related to 

20 the direction information related to the pixel. In another 
embodiment of the invention, the amount by which the 
predetermined pattern is rotated is equal to the direction 
information for that pixel. In another embodiment of the 
invention, the predetermined pattern is a split rectangle or 

25 trapezoid centered on the pixel, the rectangle or trapezoid 
having a large aspect ratio. The spiculations are identified 
by using information from the cumulative array formed by an 
accumulation of the regions of potential intersection for the 
pixels in the image. 

30 Advantageously, a given pixel in the digital 

mammogram image is considered only once in the process of 
developing the cumulative array. In this sense, a CAD system 
according to the present invention operates in the "forward 
direction, " and is very fast when compared to the "backward 

35 direction" algorithms presented previously, which use image 
pixels multiple times in deriving spiculation metrics. Thus, 
the method used in a CAD system according to the present 
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invention is very fast. Moreover, the method is highly 
amenable to hardware implementation using parallel 
processors, thus increasing the speed of the CAD system even 
further . 

5 In another embodiment of the present invention, the 

cumulative array comprises fewer pixels than the digital 
mammogram image. For example, where the digital mammogram 
image is M x N pixels, the cumulative array may be 0.25M x 
0.25N pixels. Using line information and direction 

10 information related to each digital mammogram pixel, a region 
of potential intersection is determined for each pixel with 
respect to the smaller cumulative array and proportionally 
located in a smaller 0.25M x 0.25N space. The regions of 
potential intersection are then accumulated into the 0.25M x 

15 0.25N cumulative array. Because there are addition 

operations taking place for fewer cumulative array pixels, 
the algorithm is made significantly faster without a 
significant loss in resolution. 

In another embodiment of the invention, the CAD 

20 system performs the step of computing line information and 
direction information for each image pixel in the digital 
mammogram, followed by the step of computing a weighting 
function WT(theta) based on statistical information taken 
among direction information for all image pixels. The 

25 direction information for each digital mammogram image pixel 
having coordinates (i,j) is an angle THETA ( i , j ) , and the 
weighting function WT(theta) is equal to WT (THETA ( i , j ) ) for 
that pixel . For each digital mammogram pixel , a region of 
potential intersection is determined and accumulated into the 

30 cumulative array after being weighted by WT (THETA (i , j ) ) for 
that pixel. The weighting function WT(theta) is computed by 
calculating a histogram function H(theta) of the direction 
information THETA ( i , j ) for all image pixels, followed by the 
step of developing the function WT(theta) as having an 

35 inverse relationship to the histogram function H(theta). In 
this manner, lines perpendicular to a predominant line 
direction in the digital mammogram image are emphasized, 
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whereas lines parallel to the predominant line direction in 
the digital mammogram are de- emphasized, thus increasing 
system precision and accuracy. 

In another embodiment of the present invention, 
5 mass information corresponding to the digital mammogram 
image, including mass location information, is computed in 
addition to the cumulative array. Information in the 
cumulative array is used in conjunction with the mass 
information for identifying regions of interest in the 
10 digital mammogram image, such as by using a linear classifier 
method based on mass information and cumulative array 
metrics . 

In another embodiment of the invention, local 
attention is given to the cumulative array near locations 

15 having a strong circumscribed mass candidate. The cumulative 
array is thresholded by a first threshold value in a first 
region not including the strong circumscribed mass candidate 
location, whereas the cumulative array is thresholded by a 
second value less than the first value in a second region 

20 which includes said strong circumscribed mass candidate. In 
this manner, spiculations which otherwise would have fallen 
below a threshold value in the cumulative array are detected 
when associated with a strong circumscribed mass candidate, 
for assigning a value of spiculatedness to said circumscribed 

25 mass candidate in mass detection and classification. 

In another embodiment of the invention, mass 
information for the digital mammogram image is computed by 
using information in a sphericity array. The sphericity 
array is calculated by the steps of computing a gradient 

30 plane from the digital mammogram image, the gradient plane 
having pixels, each gradient plane pixel having a gradient 
intensity value and a gradient direction value, 
determining a region of potential centroid for each gradient 
plane pixel using the gradient intensity value and gradient 

35 direction value for that pixel, and accumulating the regions 
of potential centroid to produce a sphericity array. In this 
way, strong circumscribed mass candidates may be detected by 
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using information derived from an algorithm which is a 
forward direction algorithm similar to the forward direction 
algorithm for detecting spiculations . 

In another embodiment of the present invention, the 
5 fast CAD system is capable of locating noneccentric 

spiculations for increased precision, accuracy, and reduction 
of false positives. The fast CAD system according to this 
embodiment performs a method comprising the steps of 
determining a region of potential intersection for each of a 

10 plurality of image pixels using line information and 

direction information related to that image pixel, computing 
a plurality of weights corresponding to each of the plurality 
of image pixels, accumulating for each of the plurality of 
image pixels a plurality of weights into a plurality of 

15 accumulation planes for those pixels located within the 

region of potential intersection for that image pixel, and 
processing information contained in said plurality of 
accumulation planes for identifying the noneccentric 
spiculations in the image. The plurality of accumulation 

20 planes comprises a first accumulation plane ACC1, a second 
accumulation plane ACC2, and a third accumulation plane ACC3 , 
the first, second, and third accumulation planes ACC1, ACC2 , 
and ACC3 being processed for producing a spiculation activity 
plane ACT and a spiculation eccentricity plane ECC for use in 

25 locating noneccentric spiculations. 

The spiculation activity plane ACT and the 
spiculation eccentricity plane ECC are computed using 
information in the plurality of accumulation planes such that 
the spiculation activity plane ACT comprises pixel values 

30 related to the presence of spiculations, and such that said 
spiculation eccentricity plane ECC comprises pixel values 
related to the presence of eccentric spiculations. A 
spiculation output plane SO is formed by setting, for each 
pixel (i,j), SO(i,j) equal to a first constant multiplied by 

35 ACT(i,j) added to a second constant multiplied by ECC(i,j), 
the first constant typically being a positive number and the 
second constant usually being about -0.5 times the first 
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constant. In this manner, the spiculation output plane 
SO(i,j) will contain high values near locations having 
spiculations, and will contain maxima among these high values 
corresponding to spiculations which are less eccentric and 
5 more radially symmetric. In this manner, false positives are 
reduced and accuracy and precision are increased. 

In another embodiment of the invention, mass 
information corresponding to the digital mammogram image is 
computed, the mass information including mass events, each 
10 event comprising mass centroid location, mass area, mass 

elongation, and mass contrast. Information contained in the 
SO, ACT, and ECC arrays is used in conjunction with the mass 
information in a linear classifier method for identifying 
regions of interest in the image. 

15 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. la shows an outside view of a computer aided 
diagnostic (CAD) system according to the present invention; 

FIG. lb shows a block diagram of a CAD processing 
2 0 unit of a CAD system according to the present invention. 

FIG. 2 is a flowchart representing overall steps 
practiced by the system of FIG. 1; 

FIG. 3 is a flowchart representing overall steps 
taken in the spiculation detection algorithm portion of FIG. 
25 2; 

FIG. 4 is a flowchart representing line detection 
and direction estimation steps for producing a line image and 
a direction in accordance with the present invention; 

FIG. 5 shows a portion of a binary line image 
30 corresponding to a digital mammogram image 

FIG. 6 is a flowchart representing a spiculation 
detection algorithm in accordance with one embodiment of the 
present invention; 

FIG. 7 is a conceptual diagram showing the 
35 determination of potential regions of intersection for two 
points in a line image in accordance with the present 
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invention, the two lines not being associated with a common 
spiculation; 

FIG. 8 shows a typical predetermined pattern P used 
to generate regions of potential intersection in the 
5 spiculation detection algorithm according to the present 
invention; 

FIG. 9 is a conceptual diagram showing the 
determination of potential regions of intersection for two 
points in a line image in accordance with the present 
10 invention, the two lines being associated with a common 
spiculation; 

FIG. 10 shows an image representing a portion of 
the cumulative array resulting from the performance of a 
spiculation algorithm on a digital mammogram image in 
15 accordance with the present invention. 

FIG. 11 is a flowchart representing a spiculation 
detection algorithm in accordance with a further embodiment 
of the present invention; 

FIG. 12 is a conceptual diagrams showing pixel 
20 mappings related to the steps of FIG. 11; 

FIG. 13 (a) is a simplified diagram of a typical 
line image corresponding to a digital mammogram image; 

FIGS. 13 (bl) and 13 (b2) are representative diagrams 
of eccentric and circular spiculations, respectively; 
25 FIGS. 13(c) and 13(d) show a histogram of angle 

information in a digital mammogram and a corresponding 
weighting function WT(thet a) computed from the histogram; 

FIG. 14 shows a flowchart representing spiculation 
and mass detection algorithms in accordance with a further 
30 embodiment of the present invention; 

FIG. 15 shows a flowchart representing a sphericity 
algorithm in accordance with a further embodiment of the 
present invention; 

FIG. 16 shows a flowchart representing a 
35 spiculation detection algorithm using multiple accumulation 
planes in accordance with a preferred embodiment of the 
present invention. 
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DETAILED DESCRIPTION OF THE INVENTION 

Figure la shows an outside view of a computer aided 
diagnostic (CAD) system 100 for assisting in the 
identification of spiculated lesions in mammograms according 
5 to the present invention. CAD system 100 is used as a step 
in the processing of films for mammography exams. CAD system 
100 comprises a CAD processing unit 102 and a viewing station 
104. In general, CAD processing unit 102 scans an x-ray 
mammogram into a digital mammogram image, processes the 

10 image, and outputs a highlighted digital mammogram for 
viewing at viewing station 104. 

Figure lb shows a block diagram of CAD processing 
unit 102. CAD processing unit 102 comprises a digitizer 103, 
such as a laser scanner with 50 micron resolution, for 

15 digitizing a developed x-ray mammogram 101, the x-ray 

mammogram 101 being shown in Figure la at an input to the CAD 
processing unit 102. CAD processing unit 102 generally 
includes elements necessary for performing image processing 
including parallel processing steps. In particular, CAD 

20 processing unit 102 includes elements such as a central 

control unit 105, a memory 108, a parallel processing unit 
110, and I/O unit 112. It is to be appreciated that in 
addition to the spiculation detection algorithms disclosed 
herein, processing unit 102 is capable of performing a 

25 multiplicity of other image processing algorithms such as 
mass detection and linear classification algorithms, either 
serially or in parallel with the disclosed spiculation 
detection algorithms. 

Viewing station 104 is for conveniently viewing 

30 both the x-ray mammogram 101 and the output of the CAD 

processing unit 102 on a display device 118. The display 
device 118 may be, for example, a CRT screen. The display 
device 118 typically shows a highlighted digital mammogram 
corresponding to the x-ray mammogram 101, the highlighted 

35 digital mammogram having information directing the attention 
of the radiologist to special areas which may contain 
spiculations as determined by image processing steps 
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performed by the CAD processing unit 102. In one embodiment 
of the invention, the highlighted digital mammogram will have 
black or red circles superimposed around those locations 
corresponding to spiculated lesions. 
5 Viewing station 104 also comprises a backlighting 

station 12 0 for viewing the actual x-ray mammogram 101 
itself. The radiologist is assisted by the CAD system 100 by 
viewing the display device 118, which then directs the 
attention of the radiologist to the spiculated portions of 

10 the actual x-ray mammogram 101 itself. It is to be 

appreciated that the CAD processing unit 102 is capable of 
performing other image processing algorithms on the digital 
mammogram in addition to or in parallel with the algorithm 
for detecting spiculations in accordance with the present 

15 invention. In this manner, the radiologist may be informed 
of several suspicious areas of the mammogram at once by 
viewing the display device 118, spiculations being one 
special type of the suspicious area. 

After x-ray mammogram 101 has been developed, it is 

20 inserted into the CAD system 100, which will ideally be 
located near the x-ray development area of a mammography 
clinic. After being digitized by digitizer 103, the x-ray 
mammogram will be transported using means not shown to the 
viewing station 104 for viewing by the radiologist along with 

25 the output of the display device 118 as described above. 
After the x-ray mammogram 101 has passed through the CAD 
system 100, it will be taken away and will undergo the same 
processing currently practiced in clinics. It is to be noted 
that memory 108 of CAD processing unit 102 may be used in 

30 conjunction with I/O unit 112 to generate a permanent record 
of the highlighted digital mammogram described above, and/or 
may also be used to allow non-real-time viewing of the 
highlighted digital mammogram. 

Figure 2 shows the general steps performed by CAD 

35 processing unit 102 on the x-ray mammogram. At step 202, the 
x-ray mammogram is scanned in and digitized into a digital 
mammogram. The digital mammogram may be, for example, a 3 00 0 
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x 4000 array of 12 -bit gray scale pixel values. Such a 
digital mammogram would generally correspond to a typical 8 " 
x 10" x-ray mammogram which has been digitized at a 50 micron 
spatial resolution. Because a full resolution image such as 
5 the 3000 x 4000 image described above is not necessary for 
the effectiveness of the present invention, the image may be 
locally averaged, using steps known in the art, down to a 
smaller size corresponding, for example, to a 200 micron 
spatial resolution. At such a resolution, a typical image 

10 would then be an M x N array of 12 -bit gray scale pixel 

values, with M being near 900, for example, and N being near 
12 00, for example. In general, however, either the full 
resolution image or the locally averaged image may be used as 
the original digital mammogram in accordance with the present 

15 invention. Without limiting the scope of the present 
invention, and for clarity of disclosure, the "digital 
mammogram image" is considered to be an exemplary M x N array 
of 12 -bit gray scale pixel values. 

Figure 2 shows the digital mammogram image being 

20 processed at step 204 by an overall spiculation algorithm in 
accordance with the present invention. As discussed 
previously, the overall spiculation algorithm performed at 
step 204 generates a list of locations in the digital 
mammogram image which may correspond to spiculations , along 

25 with information such as spiculation intensity for each 
location. Following step 2 04, the digital mammogram image 
and list of spiculation locations and information is sent for 
display to the viewing station 104 at step 206. 

Figure 3 shows in more detail the steps associated 

30 with the overall spiculation algorithm of step 204 of Figure 
2. In particular, the overall spiculation algorithm 
comprises the step of generating an M x N line image and an M 
x N direction image from the digital mammogram image at line 
and direction detection step 302. The M x N line image 

35 generated at step 302 comprises, for each pixel (i,j), line 
information in the form of a "1" if the pixel (i,j) has a 
line passing through it, and a "0" otherwise. The M x N 
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direction image further comprises, for each pixel (i,j), 
direction information in the form of a number THETA. If the 
line information is a "1" for the line image pixel (il,jl), 
the number THETA for the pixel (il,jl) corresponds to an 
5 estimated direction of the tangent to the line passing 
through (il, jl) . 

Figure 3 also shows a spiculation detection and 
prioritization step 304 following line and direction 
detection step 302. Spiculation detection and prioritization 

10 step 304 comprises the step of identifying, classifying, and 
prioritizing spiculations in the digital mammogram image by 
processing information in the line image and direction image 
generated by line and direction detection step 3 02. 

Figure 4 shows steps corresponding to the line and 

15 direction detection step 302 of Figure 3. Shown at Figure 4 
is the direction detection step 402 for detecting at each 
pixel (i,j) a direction corresponding to a lines, if any, 
passing through the pixel (i,j) in the digital mammogram 
image. Direction detection step 402 comprises the step of 

20 separately convolving the digital mammogram image with three 
Gabor kernels K(0) , K(60), and K(120) . The Gabor kernels are 
derived from the Gabor filter which, as known in the art, is 
the second derivative of a Gaussian kernel given by: 

25 l 

G{r,o) = ^_ ex p((-r 2 )/(2a 2 )) (1) 

2tzo 2 



The second derivative of this function with respect 
to x, quantized into a finite sized integer array, yields the 
K(0) kernel. Further, by rotating this array by 60 degrees 
and 12 0 degrees, the two other kernels K(60) and K(120) are 
obtained. The step of separately convolving the digital 
mammogram with the kernels K(0) , K(60) , and K(120) yields 
three images W(0), W(60), and W(120), respectively. 

At step 402, direction information theta(i,j) for 
each pixel (i,j) is obtained by using a formula such as that 



- 16 - 



WO 99/63478 



PCT/US98/11729 



disclosed in Koenderink and Van Doorn, "Generic Neighborhood 
Operators," IEEE Transactions on Pattern Analysis and Machine 
Intelligence, Vol. 14, No. 6 (June 1992) and given by: 



ft- l,Mn./?f fr(60) -FIT(120) ) ( 2 ) 



Figure 4 further shows line detection step 404 for 
detecting line information in the digital mammogram image. 
Positive contrast (light) lines are important, as opposed to 
negative contrast (dark) lines, since the former is how 
spiculations are manifested in x-ray films. Line detection 
step 404 comprises the step of deriving a function W(THETA) 
from the images W(0), W(60), and W(120) using a formula 
disclosed in the Koenderink reference cited supra: 

^(0) =-i (l+2cos (26) ) W o (0) +^ (1-cos (20) V3sin(20) ) W o (60) 

+ ± (1-cos (20) -v/3sin(20) ) ^(120) (3) 



25 After being computed, W(THETA) is thresholded at 

some positive value for obtaining a binary line image. 

Figure 5 shows a portion of binary line image 
corresponding to a digital mammogram image, the white pixels 
corresponding to a value "1" and the dark pixels 

30 corresponding to a value of "0". The significance of the 
binary line image is that every point in the white (line) 
areas of the binary line image belongs to a line of 
sufficient strength to exceed the threshold chosen. 

Thus, after step 4 04 of Figure 4, there exists a 

35 line image and a direction image corresponding to the digital 
mammogram image available for further processing by the 
spiculation detection and prioritization step 304 of Figure 
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3 . While there are several methods known in the art for line 
and direction detection, the above approach is employed in a 
preferred embodiment of the present invention because the 
computationally intensive parts consist of the three 
5 convolutions performed to obtain W(0), VI (60) , and VI (120) , and 
these convolutions are easily implemented in a highly 
parallel processor such as that used in processing unit 104. 
By implementing these convolutions in the spatial domain in a 
hardware parallel processor, the speed of computation easily 

10 meets normal through-put requirements for clinical practice. 
For notational simplicity, the direction information in the 
direction image corresponding to a pixel (i,j) will be 
denoted THETA ( i , j ) . 

Figure 6 shows a block diagram outlining steps for 

15 accomplishing the spiculation detection and prioritization 
step 304 of Figure 3. In particular, Figure 6 shows the 
feed-forward step 602 of generating a cumulative array C from 
line and direction information corresponding to each pixel in 
the line image. The cumulative array C, which is generally 

20 the same size as the M x N digital mammogram image, the M x N 
line image, and the M x N direction image, is first 
initialized. Then each pixel (i,j) in the digital mammogram 
image is considered. In particular, if a pixel (i,j) has 
line image information with a value "0", it is ignored. 

25 However, if the pixel (i/j) has line image information with a 
value "1", the cumulative array C is incremented by a 
constant value for each pixel located in a region 
corresponding to a predetermined pattern centered at the 
pixel location (i,j). The predetermined pattern is generally 

30 a split rectangular or trapezoidal pattern of pixels having a 
high aspect ratio, rotated at an angle corresponding to the 
direction information for the pixel (i,j). The next pixel in 
the digital mammogram image is then considered. After each 
pixel in the digital mammogram image is considered, the 

35 cumulative array C will be completely formed. 

Figure 6 further shows a step 604 for identifying 
and prioritizing spiculations based on information contained 
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in the cumulative array C. Generally, the cumulative array C 
will contain values C(i,j) corresponding to the strength of 
spiculations centered on the digital mammogram pixel located 
at ' coordinates (i,j). Generally, the cumulative array C will 
5 contain local maxima where spiculations, if any, are likely 
to be present. The cumulative array C may be processed in a 
variety of ways for identifying locations and strengths of 
spiculations in the digital mammogram image. For example, in 
the simplest procedure, step 604 may comprise the step of 

10 identifying those locations (i,j) in the cumulative array C 
for which C{i,j) is greater than a predetermined threshold 
value. However, other ways for identifying and prioritizing 
spiculations based on information contained in the cumulative 
array C, many of which will be apparent to one of ordinary 

15 skill in the art upon reading this disclosure, may be used in 
step 604. 

Figure 7 shows a conceptual diagram of the addition 
of patterns for each pixel in the line image. Shown in 
Figure 7 are two points (il,jl) and (i2,j2) which lie 

20 somewhere along lines LI and L2 , respectively. Regions of 
potential overlap for these two points, denoted P{il,jl) and 
P(i2,j2), are shown superimposed over the line image for 
clarity. The regions of potential overlap P(il,jl) and 
P(i2,j2) are simply rotated versions of a predetermined 

25 pattern P, shown in Figure 8, which have been translated to 
be centered on the points (il,jl) and (i2,j2), respectively, 
the amount of rotation being a direction image value 
THETA ( i 1 , j 1 ) for pixel (il,jl) and THETA ( i 2 , j 2 ) for pixel 
(i2,j2). As shown in Figure 8, the predetermined pattern P 

30 is of a split rectangular or trapezoidal shape having a high 
aspect ratio (i.e., a large width to height ratio) . 

As can be seen in Figure 7, because of the choice 
of a high- aspect -ratio pattern P, the region of potential 
overlap P(il,jl) for the point (il,jl) is roughly equivalent 

35 to a tangent of the line LI containing the point (il,jl), the 
tangent having been truncated at a distance corresponding to 
the length of the pattern P. The spiculation detection 
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algorithm according to the present invention is based on the 
principle that any two lines LI and L2 belonging to the same 
spiculation will have an entire series of points (il, jl) and 
(i2,j2) whose tangents will overlap near a common point at 
5 the center of the spiculation. In Figure 6, the lines 

denoted LI and L2 do not belong to a spiculation, and it can 
be seen that the regions of pattern overlap for various 
points (il,jl) and (i2,j2) along lines LI and L2 will be 
dispersed throughout many locations in the cumulative array 

10 C. In contrast, as shown in Figure 9, lines L3 , L4 , L5 , and 
L6, having pixels including exemplary pixels (i3,j3), 
(i4,j4), (i5,j5), and (i6,j6), respectively, belonging to the 
same spiculation will have repeated overlap of tangents near 
the center of the spiculation. 

15 Accordingly, if the regions of potential overlap, 

denoted generally as P(x,y), are accumulated into the 
cumulative array C, the cumulative array C will contain 
higher values at locations corresponding to spiculations . 
The greater the number of intersecting, radiating lines at a 

20 spiculation, the greater the value of the cumulative array C 
at the center of that spiculation. 

Importantly, it is noted that the spiculation 
detection algorithm according to the present invention is a 
"forward direction" algorithm. Each pixel in the line image 

25 is processed only once in generating the cumulative array. 
Furthermore, a lookup table procedure is used which directly 
maps the digital mammogram pixels (i,j) lying along lines and 
having direction values THETA ( i , j ) into regions of covered 
pixels in the cumulative array C, based on the shape of the 

30 predetermined pattern P. This provides for an especially 
fast production of the cumulative array C. 

According to another embodiment of the invention, 
the lookup table procedure for mapping points and directions 
into regions of potential intersection incorporates weighting 

35 values in the predetermined pattern P. Thus, instead of 

simply adding a "1" to the cumulative array C for every array 
location falling inside the pattern P(il,jl), a weighting 
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value may be used. For example, points corresponding to a 
point Pa in the predetermined pattern P of Figure 8, lying 
directly along the center line of P, may be assigned a 
heavier weight than a point Pb lying along the periphery of 
5 the predetermined pattern P. This is because peripheral 
points are less likely than center line points to lie 
directly at the center of the spiculation. 

The predetermined pattern P is designed in an 
empirical manner to maintain accuracy (minimizing false 

10 positives) while maximizing precision in finding spiculations 
of the appropriate size. In a prime example, it may be 
desirable to detect spiculations which have a radius of 
around 5 mm, because if the spiculations grow to a size much 
larger a 5 mm radius, it may be too late to stop the spread 

15 of the cancer. For a digital mammogram image in which 1 
pixel corresponds to 2 00 microns, or in which 2 5 pixels 
equals 5 mm, it has been found that a predetermined pattern P 
having its opposing trapezoids extending from about 15 pixels 
to about 80 pixels from the center, with a width ranging from 

20 about 4 pixels near the center to 20 pixels at the periphery, 
yields satisfactory results. It is noted, however, that the 
scope of the invention may encompass a wide variety of shapes 
and dimensions of the predetermined pattern P, as empirically 
and/or analytically determined to optimize system performance 

25 based on desired performance parameters. 

Figure 10 shows in a pictorial fashion the 
cumulative array C generated from the line image of Fig. 5. 
A side -by- side comparison of Figure 10 and Figure 5 shows 
that the cumulative array C contains a local maximum at a 

30 point corresponding to the intersection of several lines 

forming a spiculation in the line image of Figure 5. Using 
one of many methods, the cumulative array C may be used to 
generate highlighted regions for display on the display 
device 118 of Figure 1. For example, the cumulative array C 

35 may be thresholded at a certain level, and points exceeding 
the threshold value may be displayed. Alternatively, a 
selected number of top values in the cumulative array C may 
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be selected for display. Other methods may be used as well 
as will be discussed below. 

Figure 11 shows a flowchart showing a spiculation 
detection and prioritization step 1104 representing a 
5 preferred embodiment of the invention. The spiculation 

detection and prioritization step 1104 is to be performed in 
place of the spiculation detection and prioritization step 
304 of Figure 3. In the preferred embodiment, the M x N line 
and direction images have been computed according to the 

10 above described methods. However, in the preferred 

embodiment of step 1104, the cumulative array C is smaller 
than the M x N digital mammogram image, having dimensions AM 
x BN, where A < 1 and B < 1. For example, A and B may both 
be chosen as 0,25, in which case the cumulative array C has 

15 dimensions 0.25M x 0.25N. Step 1104 comprises the step 1106 
of generating the cumulative array C in a many-to-one 
mapping, followed by the step 1108 of identifying and 
prioritizing spiculations based on information in the 
cumulative array C. 

20 In step 1106, a region of potential intersection is 

determined for each pixel (i,j) in the digital mammogram 
image based on the line information and direction information 
for the pixel (i, j) . For each digital mammogram image pixel 
(i # j) having line image information with a value "1", the 

25 cumulative array C is incremented by a constant value for 
each pixel located in a region corresponding to a 
predetermined pattern P' centered at a pixel location (k,l) 
in the cumulative array C. In particular, the predetermined 
pattern P' is a proportionally reduced version of the 

30 function P shown in Figure 8 for application to the smaller 
cumulative array C, and (k,l) is chosen such that k = INT(Ai) 
and 1= INT(Bj). Thus, for example, where P extended from an 
inner radius of 15 pixels to 80 pixels, the array P' would 
extend from about 4 pixels to 20 pixels for A = B = 0.25, As 

35 before, the region of potential intersection is the 

predetermined pattern P' rotated by the amount THETA ( i , j ) 
corresponding to the digital mammogram image pixel (i,j). 
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The step 1106 is preferable to step 602 in that 
fewer computations are needed to produce the cumulative array 
C, thus saving time. This is because fewer pixels in the 
cumulative array C are incremented for a given pixel (i,j) in 
5 the digital mammogram image. The computation time is reduced 
by a factor roughly equal to the product AB, because only 
that proportion of pixels in the cumulative array C are 
incremented when compared to the algorithm of step 602. 

Figure 12 diagrammatically illustrates the "many- 

10 to-one" mapping which occurs between the M x N digital 
mammogram image and the AM x BN cumulative array, for the 
case where A = B = 0.25. The mapping is "many- to-one " 
because the same pixel (kl,ll) in the cumulative array acts 
as the center for the accumulation of up to 16 different 

15 regions of potential intersection corresponding to points 
(il,jl) . . . (il+3,jl+3) in the digital mammogram image. 

Measurements have indicated that the CAD system 100 
performing spiculation detection in accordance with the 
present invention represents a drastic increase in speed over 

20 the Karssemeijer algorithm. Indeed, it has been found that 
spiculation detection using the step 304 of Figure 6 results 
in a several-fold increase in the speed of spiculation 
detection. In particular, it has been found that where the 
Karssemeijer method takes 1-2 hours to generate an output, 

25 the spiculation detection method using the step 3 04 of Figure 
6 takes less than 20 minutes. Even further, it has been 
found that if the CAD system 100 uses the step 1104 of Figure 
11 in spiculation detection, the computation time is reduced 
to about 20 seconds. This drastic increase in speed allows 

30 the CAD system 100 to be used in everyday radiology practice 
to increase the accuracy and reliability of the spiculation 
detection process while not slowing down the radiologist. 

Other embodiments for improving results are 
described below. It will be noted that while the following 

35 disclosed principles apply to either the case where the 

cumulative array C is a full-scale M x N array or a smaller 
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AM x BN array, the case where C is a full-scale array is 
described for conciseness and clarity of disclosure. 

Figure 13 illustrates a method of computing weights 
WT(theta) for weighting regions of potential intersection 
5 before accumulation into the cumulative array C in accordance 
with another embodiment of the present invention. It has 
been found that lines in a digital mammogram of a real breast 
will tend to be directed toward the nipple of the breast, as 
shown in Figure 13 (a) . The presence of many near-parallel 

10 lines tends to cause detection of highly eccentric 

spiculations as shown in Figure 13 (bl) , resulting in false 
positives in the spiculation detection algorithm. It is more 
desirable to locate and identify the circular, quasi- 
symmetric -type spiculations as shown in Figure 13 (b2), which 

15 correspond to spiculations which are more likely to be 
suspicious, than to identify the highly eccentric 
spiculations of Figure 13 (bl) , which are often simply the 
result of lines facing near the predominant direction which 
incidentally intersect. 

2 0 It has been found that weighting the regions of 

potential intersection by an appropriate function 
WT ( THETA ( i , j ) ) for the pixel being accumulated will 
appropriately de-emphasize predominant -direction lines and 
emphasize lines near the perpendicular to the predominant 

25 direction. The function WT(theta) is computed by first 
computing statistical information, such as a histogram 
H(theta), of pixel directions THETA ( i , j ) following step 302 
of Figure 3. A typical distribution H(theta) is shown in 
Figure 13 (c) , the predominant direction indicated being near 

30 theta = 90 degrees in the example shown. The function 
WT(theta) is then determined by plotting an inverse to 
H (theta) as shown in Figure 13(d). Generally, WT (theta) will 
be a minimum where H (theta) is a maximum. 

The exact shape of the plot of WT (theta) may be 

35 varied according to empirical data for yielding the best 
results. However, WT (theta) will have a general shape as 
shown in Figure 13 (d) . The cumulative array C is incremented 
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at those locations falling within the rotated and shifted 
region P<i,j), and is incremented by an amount equal to 
WT (THETA ( i , j ) ) for the pixel (i,j) instead of a constant 
vailue, as described in previous embodiments. In this manner, 
5 pixels located along lines in the predominant direction are 
de-emphasized in generating the cumulative array C. Also, 
pixels located along lines near-perpendicular to the 
predominant direction are emphasized. This results in the 
ability to detect spiculations with greater accuracy and 

10 precision, and with fewer false positives due to incidental, 
highly eccentric spiculations. 

Figure 14 shows a flowchart illustrating steps for 
identifying regions of interest to the radiologist in an x- 
ray mammogram according to another embodiment of the present 

15 invention. Figure 14 shows a step 1402 similar to step 202 
disclosed above, a step 1404 similar to step 302 disclosed 
above, and a step 1406 similar to step 602 disclosed above. 
Performed in parallel with steps 1404 and 1406, however, is a 
step 14 08 for generating mass information. Step 1408 may be 

2 0 performed by using steps known in the art for detecting, 
classifying, and prioritizing masses. Such steps are 
disclosed, for example, in Yin et al . , "Computerized 
Detection of Masses in Digital Mammograms: Analysis of 
Bilateral Subtraction Images," Med. Phys . 18(5) (Sept/Oct 

25 1991), and in U.S. Pat. No. 5,133,020 to Giger et al, 

entitled "Automated Method and System for the Detection and 
Classification of Abnormal Lesions and Parenchymal 
Distortions in Digital Medical Images," the disclosures of 
which are hereby incorporated by reference into the present 

30 application. Mass information may include, but is not 

limited to, mass events, each event comprising, for example, 
mass centroid location, mass area, mass elongation, and mass 
contrast. Mass information may also include information 
derived from region-growing algorithms known in the art and 

35 described, for example, in Gonzalez, Digital Image Processing 
at pp. 369-375, the disclosure of which is incorporated 
herein by reference . 
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Figure 14 also shows a step 1410 for identifying 
regions of interest using the cumulative array C along with 
the above mass information for detecting regions of interest 
in the digital mammogram. It has been found that mass 
5 information may be used in conjunction with spiculation 
information for reducing certain forms of false positives 
resulting from radial lines that originate from accidental 
crossing of lines, or from objects such as "radial scars" 
from previous biopsies. Many of these false positives can be 

10 eliminated at step 1410 by requiring, for local maxima in the 
cumulative array C, or for other spiculation locations 
resulting from operations such as thresholding of the 
cumulative array C, the existence of a central mass density 
from which the lines originate. Thus, as an example of this 

15 embodiment, if there is no mass associated with a local 
maximum in the array C, as dictated by the above mass 
information, then this point will not be identified to the 
radiologist as a highlighted region even though the 
spiculation identification criteria would otherwise have been 

20 met. In this way, false positives are reduced. 

Figure 15 shows a novel alternative for step 14 08 
in determining mass information according to another 
embodiment of the present invention. It has been found that 
the spiculation detection algorithm as exemplified by steps 

25 302 and 304 may be modified for detecting masses and for 
computing mass information for use in step 1410. In 
particular, it has been found that if at step 302, the 
gradient of the digital mammogram image is computed instead 
of the line and direction images, and then the step 3 04 is 

30 executed using the gradient intensity and gradient direction 
for each pixel (i,j), instead of the binary line value and 
line direction THETA ( i , j ) for each pixel (i,j), respectively, 
the resulting cumulative array, which will now be called a 
"sphericity array," will contain high values at mass 

35 locations instead of spiculation locations. 

Shown at Figure 15 is a step 1502 for computing the 
gradient of the digital mammogram image using methods known 
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in the art, such as the Sobel matrix method. Resulting from 
the step 1502 is a gradient plane which, for each pixel 
(i,j), comprises a gradient intensity value and a gradient 
direction value. Shown at step 1504 is a step for 
5 determining, for each gradient plane pixel (i,j), a "region 
of potential centroid" using the gradient intensity value and 
the gradient direction value for the pixel (i,j). The region 
of potential centroid is computed in manner similar to the 
region of potential intersection described above, using 

10 gradient intensity instead of the binary line information and 
using gradient direction instead of line direction 
THETA (x, y) . Shown at step 1506 is the step of accumulating, 
in a manner similar to that for regions of potential 
intersection for the previous spiculation detection 

15 embodiments, the regions of potential centroids into the 

sphericity array. Finally, shown at step 1508 is the step of 
using the information in the sphericity array for determining 
mass information, using steps similar to those for using the 
cumulative array information to detect spiculations . The 

20 mass information may then be used in conjunction with the 
spiculation information at step 1410 of Figure 14, as 
discussed previously. Using this approach, the mass 
information may be computed faster than in prior art mass 
information computation algorithms by using a forward 

25 direction algorithm. 

Another embodiment of the present invention may be 
used at step 1410 for handling intermediate cases of 
spiculations. There is a class of breast masses that are not 
heavily spiculated, but rather circumscribed with small 

30 spiculations or very few spiculations which, while producing 
some signal in the cumulative array C, would not otherwise be 
detected as spiculations by the spiculation detection 
algorithm disclosed. A technique known as a "local 
attention" algorithm may be used for these intermediate cases 

35 as follows. Strong (high probability) circumscribed mass 

candidates are detected by algorithms such as those disclosed 
for step 1408 discussed previously. Once regions containing 
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strong circumscribed mass candidates are identified, the 
thresholds such as those described above for detecting high 
probability spiculation candidates are lowered in these 
regions. When this technique of "local attention" is 
5 employed in this manner, the spiculations identified by means 
of the reduced thresholds are used to further classify and 
describe a measure of "spiculatedness " for a circumscribed 
mass in order to increase its level of concern. 

Thus, in this embodiment of the invention, local 

10 attention is given to the cumulative array near locations 

having a strong circumscribed mass candidate. The cumulative 
array is thresholded by a first threshold value in a first 
region not including the strong circumscribed mass candidate 
location, whereas the cumulative array is thresholded by a 

15 second value less than the first value in a second region 
which includes said strong circumscribed mass candidate. In 
this manner, spiculations which otherwise would have fallen 
below a threshold value in the cumulative array are detected 
when associated with a strong circumscribed mass candidate, 

20 for assigning a value of spiculatedness to the circumscribed 
mass candidate. 

Figure 16 shows a method for locating noneccentric 
spiculations for increased precision, accuracy, and reduction 
of false positives according to a preferred embodiment of the 

25 invention. In general, this preferred embodiment is directed 
toward weighting the regions of potential intersection P(i,j) 
for each image pixel and accumulating the weighted regions 
into separate accumulation arrays having different 
characteristics. The separate accumulation arrays are then 

3 0 combined in a novel manner to produce a spiculation output 
array for locating noneccentric spiculations. It is to be 
noted that the method of Figure 16 as described below may be 
modified for operating in conjunction with the statistics- 
based emphasis method of Fig. 13, the smaller cumulative 

35 array method of Figs. 11 and 12, and/or the parallel mass- 
detection and processing method of Figure 14. 
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Figure 16 shows a scanning and digitizing step 1602 
similar to step 202, and a line and direction detection step 
1604 which operates in a manner similar to step 302. Figure 
16 further shows a step 1606 for generating a first 
5 accumulation plane ACC1, a second accumulation plane ACC2 , 
and a third accumulation plane ACC3 . The accumulation planes 
ACC1, ACC2, and ACC3 are each computed in a manner similar to 
the computation of the cumulative array C described with 
respect to previous embodiments, except that the amount by 

10 which the accumulation planes are incremented is equal to 
functions Wl(i,j), W2(i,j), and W3(i,j), respectively, for 
those locations in the accumulation planes falling within the 
translated and rotated predetermined pattern P(i, j) . 

Keeping in mind that the values THETA ( i , j ) were 

15 computed at step 1604, the values for Wl , W2 , and W3 are as 
follows : 

Wl(i,j)= cos 2 (THETA (i,j) ) (4) 
W2(i,j)= sin 2 (THETA (i,j) ) (5) 
W3(i,j)= 2 cos (THETA ( i , j ) ) sin (THETA (i , j ) } (6) 

20 

Figure 16 also shows a step 1608 for generating a 
spiculation activity plane ACT(i,j) and a spiculation 
eccentricity plane ECC(i,j), which are to be computed 
according to the following equations: 

25 

ACT(i,j)= Wl(i,j) + W2(i,j) (7) 
ECC{i f j) =/T(Wl(i,j) -W2(i,j)) 2 + W3 2 (i,j) ) (8) 

30 

Generally, when computed according to the above 
equations, the spiculation activity plane ACT(i,j) will be 
higher in those areas corresponding to spiculations , and both 
ACT(i,j) and the spiculation eccentricity plane ECC(i,j) will 
35 be low in those areas not corresponding to spiculations . 
Further, it has been found that, generally speaking, the 
spiculation eccentricity plane ECC(i,j) will be higher in 
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those areas corresponding to more eccentric spiculations 
exemplified by Figure 13 (bl) , but will be lower in those 
areas corresponding to more circularized spiculations 
exemplified by Figure 13 (b2) . In view of these 
5 characteristics, it has been found that a spiculation output 
array SO ( i , j ) may be computed from ACT ( i , j ) and ECC ( i , j ) such 
that SO(i,j) has desired characteristics in that SO(i,j) is 
high in those areas having circularized (i.e., noneccentric) 
spiculations and low otherwise. The step of computing 
10 SO(i,j) is shown at step 1610 of Fig. 16. 

In particular, it has been found that SO(i,j) may 
obtain the desired characteristics when computed according to 
the following equation: 

15 SO(i,j)= Rl*ACT(i,j) - R2*ECC(i,j) (9) 

In the above equation, Rl and R2 are constants 
greater than zero. It has further been found that better 
results are obtained for Rl > R2 , and in a preferred 
20 implementation, Rl = 2 while R2 = 1 . 

Figure 16 further shows a step 1612 for generating 
mass information related to the digital mammogram image. It 
is to be appreciated that the step 1612 may be performed 
serially or in parallel with the steps 1604-1610. Figure 16 
25 further shows a step 1614 for identifying regions of interest 
and a step 1616 for displaying the regions of interest. 

In general, the step 1614 of Figure 16 may be 
performed in a manner similar to the step 1410 of Figure 14. 
However, it is to be appreciated that, in light of the above 
30 descriptions, there have been described several features 
which, according to another embodiment of the present 
invention, may be taken together in one of many ways for 
indicating regions of interest responsive to the presence of 
spiculated or stellate masses. Table 1 includes a non- 
35 limiting list of some of these features. 
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TABLE 1 



Feature 


Name 


Description 


1. ECC(i,j) 


Eccentricity plane 


Inversely related 
to circularity of 
spiculations 
centered at (i,j) 


2 . ACT ( i , j ) 


Spiculation 
activity plane j 


Related to the 
presence of 
spiculations 
centered at (i,j) 


3. SO(i,j) 


Spiculation output 
plane 


Example: SO(i,j) = 
2*ACT(i,j) - 
ECC(i, j) 


4 . Sphericity (i, j ) 


Sphericity plane 


Related to 
presence of 
circumscribed 
masses centered at 
(i, j) 


5. C(i,j) 


Cumulative array 


Related to the 
presence of 

SpiLUxaLlOIlb 

centered at (i,j) 


6. Mass events 




Examples include 
mass centroid 
location, mass 
area, mass 
elongation, and 
mass contrast. 



25 A method of linear classifiers using rule-based 

cuts (thresholds) on each feature, or in combination of two 
features taken together, may be used to determine regions of 
interest. By way of non-limiting example, simply the value 
of SO(i,j) may be thresholded by a threshold value. In 

30 another example, a plot may be made of (l/ECC(i,j)) versus 
elongation for events having a value of mass area above a 
certain threshold. Minimum threshold values along the 
abscissa and ordinate may be selected, and events falling in 
the upper right quadrant may be selected as regions of 

35 interest, with a view toward not identifying large elongated 
masses, as determined by mass algorithms, unless they are 
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associated with a highly circular spiculation, as determined 
by spiculation algorithms. The values of thresholds used may 
be determined empirically by examining the distribution or 
true and false positive indications. 
5 By way of further non- limiting example, a plot of 

SO(i,j) versus (1/ECC (i , j ) ) may be made, and events falling 
within an upper right quadrant of the plot may be selected, 
with a view toward identifying large magnitude spiculations 
which are also highly circular. It can be seen that many 

10 linear classifier methods may be used based on the entries of 
Table 1 falling within the scope of the present invention. 
It is also noted that a standard back-propagation neural 
network classifier may be used in conjunction with the 
features of Table 1 and other features, the neural network 

15 classifier being trained on multiple examples of true and 
false positive indications. 

While preferred embodiments of the invention have been 
described, these descriptions are merely illustrative and are 
not intended to limit the present invention. For example, 

20 although the embodiments of the invention described above 
were in the context of a system for computer aided diagnosis 
and detection of breast carcinoma in x-ray films, those 
skilled in the art will recognize that the disclosed methods 
and structures are readily adaptable for broader 

25 applications. For example, the invention is applicable to 
many other types of CAD systems for detection of other types 
of medical abnormalities. 

30 



35 
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CLAIMS 

What is claimed is: 

1. A method of detecting spiculations in an image, the 
5 image having pixels, comprising the steps of: 

determining a region of potential intersection for each 
of a plurality of image pixels using line 
information and direction information related to 
that image pixel ; 
10 accumulating said regions of potential intersection to 

produce a cumulative array; 

processing information contained in said cumulative 
array for identifying the spiculations in the 
image . 

15 

2. The method of claim 1, wherein said step of determining 
a region of potential intersection is performed for 
substantially every image pixel. 

20 3. The method of claim 2, wherein said step of determining 
a region of potential intersection comprises the steps of: 
determining, according to the line information related 
to the image pixel, whether the image pixel is 
located along a line; and 

25 if the image pixel is located along a line, selecting a 

region centered on the image pixel corresponding to 
a predetermined pattern, said predetermined pattern 
being rotated by an amount related to the direction 
information related to the image pixel. 

30 

4. The method of claim 3, wherein said step of determining 
a region of potential intersection comprises the step of 
selecting a null region if the image pixel is not located 
along a line. 

35 

5. The method of claim 4, said predetermined pattern having 
a center, said predetermined pattern comprising two generally 
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rectangular portions symmetrically positioned around the 
center and extending from an inner radius to an outer radius . 

6. The method of claim 5, wherein the inner and outer radii 
5 of the predetermined pattern are chosen according to a radius 

of a desired spiculation size. 

7. The method of claim 6, wherein the inner radius of the 
predetermined pattern is less than the radius of the desired 

10 spiculation size, and wherein the outer radius of the 
predetermined pattern is greater than the radius of the 
desired spiculation size. 

8. The method of claim 4, the image being M by N pixels in 
15 size, said cumulative array being M by N pixels. 

9. The method of claim 4, said step of accumulating said 
regions of potential intersection comprising the step of, for 
each of said image pixels, incrementing all pixels in said 

20 cumulative array located within said region of potential 
intersection for the image pixel by a first amount. 

10. The method of claim 2, the image being M by N pixels in 
size, said cumulative array having a size of AM by BN pixels, 

25 where A < 1 and B < 1, said step of determining a region of 
potential intersection comprising the steps of: 

determining, according to the line information related 
to the image pixel, whether the image pixel is 
located along a line, the image pixel having 
30 coordinates; and 

if the image pixel is located along a line, selecting a 
region centered on a proportional center pixel 
corresponding to a predetermined pattern, said 
predetermined pattern being rotated by an amount 
35 related to the direction information related to the 

image pixel, wherein said proportional center pixel 
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has coordinates equal to said coordinates of the 
image pixel scaled by A and B. 

11. The method of claim 10, wherein said coordinates of the 
5 image pixel comprise a first coordinate and a second 

coordinate, and wherein said coordinates of said proportional 
center pixel comprise a first coordinate equal to said first 
coordinate of said image pixel multiplied by A and 
integerized, and wherein said coordinates of said 
10 proportional center pixel further comprise a second 

coordinate equal to said second coordinate of said image 
pixel multiplied by A and integerized. 

12. The method of claim 11, said step of accumulating said 
15 regions of potential intersection comprising the step of, for 

each of said image pixels, incrementing all pixels in said 
cumulative array located within said region of potential 
intersection for the image pixel by a first amount. 

20 13. The method of claim 9, wherein said first amount is a 
fixed amount independent of the direction information 
corresponding to the pixel . 

14. The method of claim 9, further comprising the steps of: 
25 computing line information and direction information 

related to substantially every image pixel; and 
computing a weighting function based on statistical 

information related to the direction information 
for substantially every image pixel and being a 
30 function of direction information; 

wherein 

said first amount is equal to the weighting function 
corresponding to the direction information related to the 
image pixel . 

35 

15. The method of claim 14, the image pixel having 
coordinates (i,j), said weighting function being a function 
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WT (theta), said direction information related to the pixel 
being an angle THETA ( i , j ) and corresponding to a tangent of a 
line passing through the image pixel, wherein said first 
amount is equal to WT (THETA ( i , j ) ) , 

5 

16. The method of claim 15, said step of computing the 
weighting function WT(theta) comprising the steps of: 

computing a histogram H(theta) of the direction 

information THETA ( i , j ) for all image pixels; 
10 computing the weighting function WT(theta) having an 

inverse relationship to the histogram H(theta), the 
weighting function WT(theta) having a minimum value 
for a value of theta corresponding to a maximum 
value of the histogram H (theta) . 

15 

17. The method of claim 1, said step of processing 
information contained in said cumulative array for 
identifying the spiculations in the image comprising the step 
of locating local maxima in said cumulative array for 

20 locating spiculations in the image. 

18. The method of claim 1, said step of processing 
information contained in said cumulative array for 
identifying the spiculations in the image comprising the step 

25 of thresholding said cumulative array for locating 
spiculations in the image. 

19. A method of detecting regions of interest in an image, 
the image having pixels, comprising the steps of: 

30 computing mass information corresponding to said image, 

including mass location information; 
determining a region of potential intersection for each 
of a plurality of image pixels using line 
information and direction information related to 
35 that image pixel; 

accumulating said regions of potential intersection to 
produce a cumulative array; and 
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using information contained in said cumulative array in 
conjunction with said mass information for 
identifying regions of interest in said image, 

5 20. The method of claim 19, said step of using information 
contained in said cumulative array comprising the steps of: 
locating a spiculation by using information in said 

cumulative array; and 
comparing said spiculation location with said mass 
10 location information for determining the presence 

of a region of interest. 

21. The method of claim 20, said step of locating a 
spiculation comprising the step of identifying a location of 

15 a local maximum in said cumulative array. 

22. The method of claim 20, said step of locating a 
spiculation comprising the step of thresholding said 
cumulative array for identifying spiculations near locations 

20 where said cumulative array exceeds a threshold value. 

23. The method of claim 22, said threshold value being a 
constant value. 

25 24. The method of claim 19, said step of using information 
contained in said cumulative array comprising the steps of: 
identifying a location of a strong circumscribed mass 

candidate using said mass information; 
thresholding a first region of said cumulative array 
30 with a first threshold, said first region not 

including said strong circumscribed mass candidate 
location; 

thresholding a second region of said cumulative array 
with a second threshold, said second region 
3 5 including said strong circumscribed mass candidate 

location, said second threshold being lower than 
said first threshold; and 
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if said second threshold is exceeded at said strong 

circumscribed mass candidate location, indicating a 
region of interest at said strong circumscribed 
mass candidate location. 

5 

25. The method of claim 19, said step of computing mass 
information comprising the steps of: 

computing a gradient plane from said image, said 

gradient plane having pixels, each gradient plane 
10 pixel having a gradient intensity value and a 

gradient direction value; 
determining a region of potential centroid for each 

gradient plane pixel using the gradient intensity 
value and gradient direction value for that pixel; 
15 accumulating said regions of potential centroid to 

produce a sphericity array; and 
using information contained in said sphericity array for 
determining said mass information. 

20 26. The method of claim 25, wherein said step of determining 
a region of potential centroid comprises the step of 
selecting a region centered on the gradient plane pixel 
corresponding to a predetermined pattern, said predetermined 
pattern being rotated by an amount related to the gradient 

25 direction value for the gradient plane pixel. 

27. The method of claim 26, said predetermined pattern 
having a center, said predetermined pattern comprising two 
generally rectangular portions symmetrically positioned 

3 0 around the center and extending from an inner radius to an 
outer radius . 

28. The method of claim 27, wherein the inner and outer 
radii of the predetermined pattern are chosen according to a 

35 radius of a desired mass size. 
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29. The method of claim 28, said step of accumulating said 
regions of potential centroid comprising the step of, for 
each of said gradient plane pixels, incrementing all pixels 
in said sphericity array located within said region of 

5 potential centroid by an amount corresponding to the gradient 
intensity value for that gradient plane pixel. 

30. The method of claim 25, said step of using information 
contained in said cumulative array comprising the steps of: 

10 identifying a location of a strong circumscribed mass 

candidate using said mass information; 
thresholding a first region of said cumulative array 
with a first threshold, said first region not 
including said strong circumscribed mass candidate 

15 location; 

thresholding a second region of said cumulative array 
with a second threshold, said second region 
including said strong circumscribed mass candidate 
location, said second threshold being lower than 

20 said first threshold; and 

if said second threshold is exceeded near said strong 

circumscribed mass candidate location, indicating a 
region of interest near said strong circumscribed 
mass candidate location. 

25 

31- The method of claim 19, said mass information including 
mass events, each event comprising mass centroid location, 
mass area, mass elongation, and mass contrast, said step of 
using information contained in said cumulative array 
30 comprising the step of using linear classifiers for 
prioritizing said mass events. 

32. A method of detecting spiculations in an image, said 
method being capable of locating noneccentric spiculations, 
35 the image having pixels, said method comprising the steps of: 
determining a region of potential intersection for each 
of a plurality of image pixels using line 
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information and direction information related to 
that image pixel; 
computing a plurality of weights corresponding to each 
of said plurality of image pixels; 
5 accumulating, for each of said plurality of image 

pixels, said plurality of weights into a plurality 
of accumulation planes for those pixels located 
within said region of potential intersection for 
that image pixel; and 
10 processing information contained in said plurality of 

accumulation planes for identifying the 
noneccentric spiculations in the image. 

33. The method of claim 32, wherein said plurality of image 
15 pixels comprises substantially every pixel in said image. 

34. The method of claim 33, wherein said step of determining 
a region of potential intersection comprises the steps of: 

determining, according to the line information related 
20 to the image pixel, whether the image pixel is 

located along a line; and 
if the image pixel is located along a line, selecting a 
region centered on the image pixel corresponding to 
a predetermined pattern, said predetermined pattern 
25 being rotated by an amount related to the direction 

information related to the image pixel, 

35. The method of claim 34, wherein said step of determining 
a region of potential intersection further comprises the step 

30 of selecting a null region if the image pixel is not located 
along a line. 

36. The method of claim 35, said plurality of accumulation 
planes comprising a first accumulation plane ACC1, a second 

35 accumulation plane ACC2 , and a third accumulation plane ACC3 , 
wherein said first, second, and third accumulation planes 
ACC1, ACC2, and ACC3 are capable of being processed for 
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producing a spiculation activity plane ACT and a spiculation 
eccentricity plane ECC for use in locating noneccentric 
spiculations . 

5 37. The method of claim 36, said plurality of accumulation 
planes, said spiculation activity plane ACT and said 
spiculation eccentricity plane ECC each comprising pixels 
with coordinates (i,j), said step of processing information 
in said plurality of accumulation planes further comprising 
10 the steps of : 

computing said spiculation activity plane ACT and said 
spiculation eccentricity plane ECC using 
information in said plurality of accumulation 
planes such that said spiculation activity plane 
15 ACT comprises pixel values related to the presence 

of spiculations, and such that said spiculation 
eccentricity plane ECC comprises pixel values 
related to the presence of eccentric spiculations; 
forming a spiculation output plane SO, said spiculation 
20 output plane comprising pixels with coordinates 

(i,j), said spiculation output plane SO being a 
function of said spiculation activity plane ACT and 
said spiculation eccentricity plane ECC at each 
pixel ; and 

25 using information in said spiculation output plane SO 

for identifying the noneccentric spiculations in 
the image. 

38. The method of claim 37, wherein SO(i,j) is set equal to 
30 a first constant multiplied by ACT(i,j) added to a second 

constant multiplied by ECC(i,j). 

39. The method of claim 38, wherein said first constant is 
equal to 2 and said second constant is equal to -1. 

35 

40. The method of claim 39, the image pixels having 
coordinates (i,j), said direction information being an angle 
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THETA ( i / j ) , said plurality of weights corresponding to each 
image pixel comprising a first weight Wl(i,j) for 
accumulating into said first accumulation plane ACC1, a 
second weight W2(i,j) for accumulating into said second 
5 accumulation plane ACC2 , and a third weight W3(i,j) for 
accumulating into said third accumulation plane ACC3 , said 
first weight Wl(i,j) being proportional to 
(COS (THETA (i , j ) ) **2 . 

10 41. The method of claim 40, said second weight W2(i,j) being 
proportional to (sin (THETA (i , j )) **2 , and said third weight 
W3(i,j) being proportional to 
2 cos ( THETA ( i , j ) ) sin (THETA (i , j ) ) . 

15 42. The method of claim 41, said step of computing said 
spiculation activity plane ACT and said spiculation 
eccentricity plane ECC comprising the steps of: 

for each pixel at location (i,j), setting ACT(i,j) equal 
to Wl(i,j) + W2(i,j); and 
20 for each pixel at location (i,j), setting ECC(i,j) equal 

to SQRT ( (Wl(i, j) -W2 (i, j) ) **2 + (W3 ( i , j ) ) **2 ). 

43. The method of claim 42, said method further comprising 
the steps of: 

25 computing mass information corresponding to said image, 

said mass information including mass events, each 
event comprising mass centroid location, mass area, 
mass elongation, and mass contrast; and 
using information contained in said SO, ACT and ECC 

30 arrays in conjunction with said mass information 

for identifying regions of interest in said image, 
including the step of using linear classifiers for 
prioritizing said mass events. 

35 44. An apparatus for detecting spiculated lesions in a 
digital mammogram, said apparatus comprising: 
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a first memory for storing an image comprising pixels, 
each pixel containing line information and 
direction information derived from said digital 
mammogram; 

5 a processor for determining a region of potential 

intersection for each of a plurality of image pixel 
using said line information and direction 
information; and 
an accumulator for accumulating said regions of 
10 potential intersection to produce a cumulative 

array ; 

wherein the spiculated lesions in the digital mammogram 
may be detected by processing information contained in said 
cumulative array. 

15 

45. A computer- readable medium which can be used for 
directing an apparatus to detect spiculations in an image, 
comprising: 

means for directing a first memory of said apparatus to 
20 store an image comprising pixels, each pixel 

containing line information and direction 
information derived from said digital mammogram; 
means for directing a processor of said apparatus to 
determine a region of potential intersection for 
25 each image pixel using said line information and 

direction information; and 
means for directing an accumulator of said apparatus for 
accumulating said regions of potential intersection 
to produce a cumulative array; 
30 wherein the spiculated lesions in the digital mammogram 

may be detected by processing information contained in said 
cumulative array. 

46. A method of detecting spiculations in an image, the 
35 image having pixels, comprising the steps of: 

computing line information and direction information for 
each image pixel, said line information being equal 
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to a nonzero constant if said pixel lies along a 
line, said line information being equal to zero 
otherwise, said direction information being equal 
to the angle THETA of a tangent to said line if 
5 said pixel lies along said line; 

determining a region of potential intersection for each 
image pixel using line information and direction 
information related to that image pixel, comprising 
the steps of: 

10 determining, according to the line information 

related to the image pixel, whether the 
image pixel is located along a line; and 
if the image pixel is located along a line, 

selecting a region centered on the image 

15 pixel corresponding to a predetermined 

pattern, said predetermined pattern being 
rotated by the angle THETA for that image 
pixel ; and 

selecting a null region if the image pixel is 
20 not located along a line; 

accumulating said regions of potential intersection to 

produce a cumulative array; 
processing information contained in said cumulative 
array for identifying the spiculations in the 
25 image, comprising the step of identifying local 

maxima in said cumulative array; 
wherein said predetermined pattern has a center, and 
wherein said predetermined pattern comprises two generally 
rectangular portions symmetrically positioned around the 
30 center and extending from an inner radius to an outer radius, 
said inner radius and said outer radius being determined 
according to a desired spiculation size. 



35 
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msm 13 ^9 btPi 8»g-r^# nut* 2 oowt 

Wcfctimmfcm t ±3Br¥W6v&a£t izm 
5A^7-u ^rnktz ztizxK) ttmftcox^ 7 v 
mx-h^x, 

--9 1 ' «s s- * co^mmmco 2 #<rc 1 ot^rs 

u 

fc£o* wc , «S8 Wtcsaw 4 T^-tfESM&WJfje 
4 sy^fcfjfrf s^y^7xnmi^ 
msm 2 ] ±Mi<?>=mi' -^mmco^yy u y ^ 

l mTfcrrs i a Ui£ k *#fi&k-t & > ft*Jf l IB 

cii^ji 3 ] ±m<7)=fmm&<o^yyv y?% » =f 
#1^— rms^i?-5<?-y y^k RiwuTpa i^:L 
i k k-r & » mm 1 lemo^ - u xsmi/m 

■? h ^rWV -fKffiS&fiWz k & & s&BOfWh k £ £ jS 
kL£<rk£=ft®kiM>. ft#JS 1 7 - U xSQfc 

«rt tt * ^F* p — 9 s 3£t ; — * ^bwjb^ a y jK u 
A - h^T»ffi-ri»j;a £ Lfci k ZWtikt-th . If 



[ 11*116] JJB<0J^0ffifflttJB&* ja)«Wfc:S3» 
^•^'n^aX^kUcik&#Sk-ri». If^iIUBW 

[|f*JM7] JJB<01flS(Qftih0Bf*t. T*y-^ 
x-^<7)l «OT4Si<K<so«#te.J: Off d i a KLfcifc 
*«f«i:-r*, IS^6fiE»^7-yxS3ft4HK«S^-r 

[ If 8 3 ±IBcOFJtS(OPJI h nfc IWBHHW3 fc: t5 ft 

Cl^v ^^-if^^ilBHffS^fiffltt®^* 
4 k ft-S ^ 4 5 y h X 0 iz Ltz Z k £ 

mmt-tz , m *r 1 e»*> 7 - y x^^^o-r- 

[ff*3S9 ] ±fBOFJt^Pil^ix/-cBtPBlfI«rtf-*3ft 

KUJtik*1#»i:-r4. lf*Tllffi«^7-yx^m 

C If 1 0 ] ±E<7)F/f5gO|® &ft^ii$IRMttrt&& 
Jt 4 y — 9 i 3K'T — ^ B$HI«HV> x U -y b 

— ¥m&jfi±mftfe<v\mim\ l zt£& tmfe%ti& ? 
a ^y^%\m-hi.ot,zLtzzt^mm.t-tt, mm 
11 1 imco 7 - y x^m^fiw x- ^ jsaa^rs . 
[it^iii 1] ^wv-^%^mw±im^\m^ 
mk%&kW£%ti&?J syyiznfe-t&J y?y x 

£dlzLtzZt*1&mt-t&, lf*«llB«^7-yx 

T^P— 9****tS-fr. ilBT^y— rMW)3ftB@E 
ft k JJBT *3tO)?SJK k t=36o*v »t^y^7xnm 
kti 0±fEM^x^b;w*#^> iati-jtv^ 

TT^V- 9*367*-^ k LT»«^S-^T, IST^y- 

f*f'-^«ty7'J y^"^4 5 y?\,zlfimth94 5 
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y?'x=FW-im& * mmv- >y° y yyt. zti$> 
izm^x » Hjn«tessw-* = F*i'-- -rm&rajg 

tfOffifflttJRfcrS: S k ft & * -f 5 y ^ £3tffi L fc± 

±IB#^c7)3eB&H k JJB-f y 9 7 x c?^ 9 Affik (c*cf 
vvt, ±fB<^^7xn^'5ASro<&-'f ^7in 

[0001] 

[00 02] 

y * 7 x a 7 '9 -A £ 7 - U xg$ftf 4 i k tc X 0 * <7)X 
b «t 3 (c Ltzy - y xgajMftfcSfciflS* i 

9 jffl *vc v * & # . 7— y xSd&SBBStefc v >T 

-AX7°y -y^fci:D2^03EmteJHt£>*U iftfe 

•nut. ytn^commt Lxmm^titz^ms.s.^ 
/7A&7-'j x^m-rtL«" , atESrv * ummcomw 

k Lxm%Zti?zffim-t : &hlb±§m3£cr>x'<.7 h)V 

[0003] m-t&iz. 7-yx5Sfcfl-3te»fct3\,vt 

*/. ax = l/<2 Ad) 

4J^^«*sis£SR v max <•*-*£«: . Mm 

¥-y^&dZ;hZ<LZclftH£%t>%:\\ *»< L/C, « 
3fc»2&t:«yf-A*^/hS<'t63tftt=. «e*j;y 
fflV»£*lT^S H e -N e U—fJtX 0 tjfcROM^l' 

7 - y x^m^fsT-coffiffl teas t^*£*^)ffiv 



( ^muraii ) xm o as L^amsftSEfcaiE lx a v 
[ o o o 4 ] tx , Wh<r)v - y x«frMrr«i. 

tfaL He-NeW1f3KttfflE<0t~AX7 , Uy^» 
&V^4j3!l^tWvXTy 7^ti0 f3ififc7> 
it. Ztihcou- i^^Sr^tim. |frlBc^7tS-&Mk|BI— 
co^SS^Sr^-T 4=S-«iJW3tSS* ji3S$-fr^±T'SV ^fc=F 

lx^&. -f^h-h. Twns—VKMmimmizML 
xjEwmmz^ikr h <dx\ m=f^—ftms.ti t m 

fiSI*** Bfev ^^mtct5 v ^X ft coffin <5> jEtOfitcflJ 0 
ffcS±#-€n?nXjii (sin0fci3VvC. 6ifi2nnX 

[0005] 

P—1f3t**lffl LT^tii-r4 J; a tzLizft&coy- 'Jx 

[ooo6i§fet, mmmz^wrth^wv-^^m 

b>^tel^-if3t^K£A^T&5£$*U Kl^-ir* 

7- y xmmK^mymmmmnmv >y -f- 

i4. 3eSHf <y^-A 8 izmtx&coft l t:J; ->t^S 
[Rl] 

5^1 

[0007] &?z. ±#-fe*n^nx^Sn^.-C§^>(cT 
«® te*J v ^ t Icoffl^ ft OfSt W 0 ^> j5C ) X h =F 

t&BNS&ytyy'v y^fefrv\ -r^^-feT^v-if* 

<0StS^> 1 / 2 cOfflffi-CT^tSSE^ail^ t , 3tKHt 
>y ^-Sr^co*!^ 1/2 ft:OM L X M.^M^MM % ~X 
^<-ttk^-otzn^ti^t^fih. LfrL&&$>. Z\ 
aims. 9\&f-*9V >yW±cDtzMzL% 
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VMIi; LTbXxU ^X«J££I&5£LT^&PifiU:. ± 

M£0:M Sy/t^fW*t, ±#H£u?nx 

*§-£-c7) J: o ZcsEWQ 9 -4 5 >9"TMfeZft5 c\ttfX'% 

[00081 K LX. mtH> #IB¥2-27 2 26 
^jm. - 3 1 9 2-t&HL WWT7- 3 5 6 

8 8#&*L %iWlW-7 - 2 8 6 9 0 2#^fg#^li> W 
tgcoJ: dSrlSHHSrJBSfeLJ: a h Dta^«7-'Jxgi 

Lxmmizy-v ^mzm vx 2-^x^9 h/vzM 

ttJU c\ti^cr>2^cox^9 h/Uco^-^th^t^zX 
KiMftkW&X^? XolzLfzy—V^Mft 

■7*— ^ MH^Tii . ±#Mi' a ? nxj&cWTl&K § fi 
tz'SM t J: : f~9¥\ (^f^7xnm) 
^fttfsx-^? YsVt. Tl^b'u?nXj5«^Tl&&£*i 

hjut zmfr-sfoitxmmtt^? 

ZmtXokzLX^h* ZcoXolzLX, %m.cr,AV9 
7xn/7Af-? * tz S V ^X-^. ? h 
/l^tSiatcLTV^. L^L&t^. rco-r-^® 

S05n*TtA, PI-BtPBll*!^#^tll.-f V?7 *U?"7A 

T^RfctE^*^ 2 as x "3 *>^< -rs i i: {i-es 

fcT^n^nXj&cDW^tM y^7xn/5Af-^ 
Jfyr'jy/ts^t, *M7'J 77707/11! 

®"pr& t b h\,z^cr)axYifim<^ht\^tzmmW 
[0009] -1^#IB¥2 - 2 7 2 2 6 -ti^fi^te . § 

n^nxjSSt'f y?7xn/7Af-^^fy7yy 
:/U R-B#p H WC?£*0 3f§ (fti^linf) CD-f ^ 
^ 7 x n^5Af-^S»S J: die Lfc^-^ffia^ 

>fy^7xn/7Af-3'&2') (iSl^n-l 
ffl) fclCii^'So (JbSVMinffl) fO-r-^W (4 



yj»7x0^7A) 0 s #7-'-^?iJ (^fy^7x 

vMinffl) c?)x^7 h;PS:-o< 0. ;ti£>c7>x^? WW 
cOT^&fc DftiHW&X^.? YlV&'fthXo 

[0010] U&>US:j&<$>. JKO-r-^HSlSrsWi, 

*CO%^<03tSIII t° »/Wl/3 ( #> S 1 / n ) <T> 
EMKM y^7x P/5Af-^ Xodz LX 

V^BBflLL, im<D\/4M.b, 3«(*>*Winffl) 

A77'j7^h 3H (&4V*lin«) OWiFBttWS 

[00 l 1 ] 3 6«c, dcOx-^^MTf^Tli. *mni 
CMteBUt 7f«0 1 /ncoralHTM y^7x0^5Af 

t-* o . *<o n mcof—? * o r ms^-c-^fiS; t^-fmtf 

[0012] 4^, mrm-1- 3 5 6 8 8-9-^82Wm 
WPF 7-28690 2^m^Zh . WJHWtcSaH-ft^F 

f-^wyr u y^ammfi 1 a i a t t^7-y x 
^\ -rti^t;t. **Wti4. l>riEc 7 )WM¥2-2 7 2 

2 6^|g^Pi^§n^^-^»l77^*«-&i:l3]«cr) 
[0013] #^T6 - 3 1 9 2 ^4&Wtett, P-# - 3fc 

7Af-^^ty7"'J y^t. *fc*5ISSafc*l*l±S-a: 

C07— VX.^Mft%$sl>Z&\^Xli. A y?7xD/5A 

v^tcoeo, dix^cD-f y^7xn/5Af- flzML 
x-^txy-u^m^ffdXd^zLx^^. ttfc 

•h, A y^7xn^5AT ; '-^Sr&nji*3*tnffl«x 
- ^ ?ij 7-U x^m ^t^S, i cn>tz 

\^cox\ ^mmzn^tLZ^y-x^? hMz^—x 

[ 0 0 1 4 ] $ ^Z, ZCDim^ 6-319 

imi*titzy-V3&mftft&x*ix. <x<7)X o^mmh 
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«lf&o$. jniltKtW&flWfiil^lRais was 

[0015] *wn&. im&w&mzm&t&fcto 

[00 16] 

n 6 i^-if3ts: **i<e*u ±ia5sa&n t m-<rmmm 

( a ) nwk^<^^Mmmco2^ 1 £TF 
^mlx^^-^?-? t Lxwm-r m=F 

*K U ( b ) Wifci«8HB«rt ft* 

i— i*%rf-? izm^x . mmmz^m-h^-mv- 
^m^mmMcoimwz^i, bMM%ti& ?a $ y 
ttz±x\ ±nmfflw%i<zttm-?mmnffl% 

bbi>lz. ±IH-f ^7xo/9Af-^fc 

m^x±m?4 sywzttm-t&j y^7xn/7A 

ffi£*#>, ( c ) ±I3t$5£«03fcRi£i:±i5-f ^7xO 
[0017] Z<ny-VX$&ft%M(07 t --Sr!mjjm 

<e*u mmmz^m-t^T»^~^mmco^mmm^> 

1/2 J; 0 i$SV^ (fltf.tr* 1/20) IfflBittV 
Tl> y/$tu ^i??>vmk%tiX 3 r&V-<fytr : -? 
Rlf^mtT-? b LXWmZith . 0 . T^l^- 



[0018] -etT. ic?)J;dtcSW§^T*^-if' 
b TifcSfcr- 9b*. flRtr^ytfA-^Srffl 

^h-%b° it, v 7 h fx rwtc»ri> i^t-io. 

^y^7xn/7A#$»i5^. $r*5, flRf±\ WBH 
T2-2722 6#&fR, ftWPF 6-319 2^fg, 
«fBPF7 - 3 5 6 8 8-9£«** V*te#BH¥7 - 2 8 6 

9 o 2^mi l zm^iix^^m^y-v x.^m^ 
mfemz-% ~otzc\b*maL. z. com&x'^mxsmz 

Mfet& £ tc <t 0 A y? y xn^Af- ^ Sr**. 

h ?-?%kmumxn , m^.<omm^mx'^co=F»u 
~^^-?Rv^&^~?*m&.LxmL. z.fi 
zm&mz V7 N7x TWtc^a-r-i. - 1 1 

7xO/7Af-? tCS-CS'V ^^y^7xn^7A^# 

[0019] J: 0^*Wtc«. 4-f » FJrSo|3gA,tL7tB# 

(Bt^j) wwbhtih, z.zx\ wmmz^mth^w 
u-?'%3m b Mm b immmmzfo &cox\ ±sm 
fflvmizttm-t&mizcDyemm. mmi-fi&±Mz?J s 
y^zn^-t^^^mmtm^z%^fih. zti 
bmmz, wm § ntz ?izm^x, ±is^ 

fo&<7)x\ z\<mmmmzgr3\^x±si.94 sy/ (w 

7xD/7Af-?MffKk ^^xA»o^y^7 

^ tC* "CfV ^^^7xn /7AW < 

[0020] *»fflt:3&>*»4 7- ya:3»4Hfc»<sof- 
if s ^*C0*^^SSMt° 'y^CO 1/2 cot" >y -f-TM >- 
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JJEJ5r}£<0G»lttJB f 3IS5lflW>JJH*n ? a 
* SRJR-r S i k £ 4 . § fcf i . ±f Bm^«f5ffl« 
t7f^y^7in/7Af^ £«Kf- 4 £ 2: AT 

3*l4<o^ M7-yx^3B^^M»JK*^ Lte 

[0021] _kia 7 - u x^^^^^-f- ^Ma^-ffi 

ism t % h j: -5 *= a/d smss^aw^ptt k =vm<rffi 
mmmmt itmmh £ k . 

[0022] JJB7- yx^^3t^^x-^M 
3S*rr 4 =f& V —^f^mSUfiWz k & h &BOSWh k & 

Ua-bSJ* Saviztky-Golayfe) Titffi-r& 

[00243 ±ia y-v xsaawHisio-r- **QH:£r8; 

[0025] JtJE^-yxS^HfcSO'r-?*!* 



^y^mismmmrgitffrfr&if*. z-coxoizlx i 
2<kfo&^3&%t*<v%m^&VAZftii$y < v 

[0026] S ^c. ±157- 'J x^^^Ot-"-^ 
4^-f Sy^^ffi-rSiot-LTtJ:^. i<0*^. 

mmm^snbhb^-otz^m^m^xhx^, =f 

[0027] =3rfe, ±K7-Ux3fi»4Hfc»(^-^«l 
H^-ft (C4J V « , T^l— If' *5SS* i "JblB^^'(4ffi 
«5JBk**k«S$*iS^'f Sy^fcMJ&rS'f 
xn^5A«, M^-f $>7tti£Xy->y°V>7Zti 
tzA >?y x Af-? (cfflnttSt&lt tt ICffi^ 

[0028] $ fefc. *^HJ^2^«««, b 

-3^X4 y?y xa^9ASro< 0. IS-f ^7in/ 
9ASr7-yx^-r4ikfcJ: y±fBH3t^x^^ h 
i a fc^r-o t ^4 y - 'J x2a&W«<DT r - # 

mwmmzti^x. (a> T^-ifJBSMEt-f-^aB 
^«o2^<o i OT^WfSsrat-mni^^ry 

Cltt4>Sr-rv?^^t-CT^^-if3K-r-?kt 

-f y7- y y itife ^f^7;Hl u^y^7x 

n/5Af-? k LTWt-TS-r — Kk > 

( b ) mfenm&niz&^mmzmt&T&u-'fx. 
izm^x , mwmz$m-t&=F®u-v s %mm 
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htbUz, _LfB-f y 9 y x n/7Af-^Kl^i ^X 
& ^- i: . ( C ) ±f B^^S^ t ±IB4" 

[0029] .r<7)7-Ux^^f|«'r'-^®Jl^a 

# SWT, 17- U x£W3tf§^«A#iJ^R£AfS 

<DX\ 117- U x^^^cO«W»St-^fc% S / N 
Jt£Alf fc:i»J_L$-fc!r& it^'l * . 
[0030] 

S.Hl iZTpcti 3 fc, r»7-'J x«fr*#t£*> v * 

2*aau3ta. is iu>x3 izx ^m^tcx^-yy 

jU41izmt%ti& . ^ryy°)VA-&m^Ltzmt. miK 

mzttWlZtiZ. -*-<^2Ff8i:S«W»9tci*i*» 
v\ flfcfrtftfBKliilJgit 1 0 fciifjfrd . ZtiZ'tim 9 . 
1 0 "CRItSixfe^jBEtt. tr-AXT U y * 8XSS 

yXl lTm^tlfzm, SS2RI*$5-1 2X'MIf$ 

[ 0 0 3 1 ] m 1 £3tS 5- 5 cr>*#fflt,Z&K 1 4 #H 
ftt>tlX\i&. ^LX. He-Nel/-ffgl5^ 

±Mtt 1 4 znmmtt>nfzm. mmm 1 frh<r>ym. 
(Dim*&mzm&£oizLxm2]/yX7 zm-oxt: 

-AXTU >y ^ 8 tAO s - Z.X2^>cr>U—*fmz-!ri-m 

flit 9 s 1 0*CRSt§*Ufc#l'- fit b'-AXT'J 

icOT^^-^tt, M3UyXl lSrM5i§LTlS2^ 
Sf59-1 2fcm*»3. jg2^it$5>-l 2»WK 
tt^ 1 7 j^BHtt ^iXTfc •? , ±m=f»\y-^mi^<^j\ 
1 7£jEiLXl/— if^-fe^lM 8 (i/tj3y7t 



mm^- (i^-ifT^ft-f-) saw-*-*. 

[0032] f£iMI9i4, »»»«9tA*3BKfc«er 

jtflBiP- iTF^^jaaHWtJ J:* lOkHzfc* 
&£5t,zft&L<®ffl%tiX^&. ztx. mm9iz 

IW»9 W^fflfctts ±fBm3£«^5£<D**fig£iEl 
Vtzc\t tefflsithy* YAy? yy°?WmihtiX\^ 
& . ±f B**PKi«f 4"? ^^o^-^tcJ;^ 9 ct)^ 

[ 0 0 3 3 ] H2«, Hi lz^y-Vx&mft%M/& 

yfy^xiyy&m^ttf, ztizix, mwm^&F?- 
m: m&x-y- ox.$mt&Hii!rizm&-t&mmm 
mztttx) mim^2obm2mm^2 ltizx-o 

Xm^fltzm. 2f- J r^*;^«A/D^mtl2 2fc 
X-ftZtlh* i^A/D^m#ff2 2^{4BUlB^7 

xhJyfJTffr^ftJiZti&f't s^fi^-x^A 
5 yfimxzfflffim^t LXA/D^m^WMct 

X'foK), Mitf BURR-BROWN|ll|(7)D S P 1 0 
2^*»fflV^il&. ££T, A/D»2 2<iDf 

(CH. A) fcttl'— 1f=F«Nl*6*A*SfL. flfi 
Jf^y^B (CH. B) ^{4-^ ^7xn/5A| 

[0034] A/DS«2 2J4s y^^m 
2 6l,z£^Xm.feZtlZ>?4 Sym-jT2 0 0kH 

z ery-yy u yyj$fflx\ ^y^ARx/^ y^v 

B^A*SiUt7-*-^^ 1 8b' 7 bA/D^mtT. D 
S P 2 3 ( f-^y^Wot : Mx(fT I ttM 
TMS3 2 0C3 1 ) fcMS. Lt, ;tl4«fy^ 
^ ( 0*T#fiW- W^^7xn^7 
Afl^'-^ ) !±DSP23 (CJ: 0— B^fWfc^^U 2 4 
fcfl»ftS*L*. *»< t-X. #»M9(01HItfO^X-2 5 
0 0 0 ^Wf-^^'A/DSi LT ^ ^ 'J 2 4 
*U>. i«f^ DSP23(1 *^)2 4l,zWfc&KX 

mntx . %nmmm$:±.iiL3 ytx-? 2 7 1 co^ 

[0035] JilT, D S P 2 3 ^ioTHfT^ixS'r- 

ft&mcommx^z^ ^xmt< mmtz . zcoy-v 
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2 0 0 k H zT"^V7°y y^'$^T^-g» . Ufo^T , 3 >*'l> a- h«J|:££fft£ £ t i Off £> 

m-^mm mmmi£mmk : &wx'fc&. 

f k =b 2k X2 + b lk X + b 0k 5^2 

[0036] ZZlZ&^X. MZ-im3lZ7K-t£d%U &U-1-'ffim$:$iLX\^. ThSr/b*>. H3*T{4. 

-^sm^-^mz^xte. 5^2 ffob 2k . b lk mmRvmmxw&^tfhtiZ'&^-i'tf ( x±. y* ) 

b 2k = (S 2 F 0k -S 0 F 2k )/ / (S 2 2-S 0 S 4 ) ^3 

b lk = F lk /S 2 ^4 

b 0 k = (S 2 F 2k -S 4 F ok )/(S 2 2-S 0 S 4 ) 

+m +m 
Sj=Sl' F, l =2i l y t H-j 
i=-m i=-m 

[00 37] *>< LT. kcOffiSrP— iFT^fi-f-x-? ^ ( x k ' , y k ' ) #'#£>fil> . x k 'W:b°-? 

( 1+m) fflB*^ (2 5 00 0-m) HS^T^ ^XJ4# h i^^Bt^J (*-f SV^) SrSiU y k 'J4b° 

(x k ',y k ') = (x k -b lk /(2b 2k ),b 0k -b lk 2 /<4b 2k )) 5^6 

fc#£U5rVtffr£fc#ifc£. Mx.fi\ 04a. x k 'Offi T'&O. ftT**U* e-^jSTft*. £*U=«fc9. 1/- 
#±IEIHJ*rtfc#^£^<Z>— H5i4. x if' = F^«^fc' 1 -^^ (b*-7ffi) (sKhA 

#<9*£- H5) fi, b-^X(4tf>Aj£i4 hA^ttt^-fSStC^r&^-C. b-^fctfbAj&fc 

v*. zzx\ mf-^ii, w&cvtmximmb iznhtix v^*^>£^x 7?ti.;t #-es & . - 
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Abstract of JP 5026919 (A) 

PURPOSE:To obtain a spectrum analyzer of which 
the measuring level accuracy is improved by 
correcting a scallop loss being inherent in fast 
Fourier transform. CONSTITUTION:A spectrum 
analyzer is equipped with an analog-digital converter 
9 subjecting an input signal to analog-digital 
conversion and a memory 1 1 storing temporarily 
digital data obtained from the analog-digital 
conversion.; Moreover, it is equipped with a 
multiplier which multiplies the digital data stored in 
the aforesaid memory 1 1 by a Gaussian function as 
a windowfunction, a fast Fourier transform unit which 
applies a fast Fourier transform processing to the 
result of multiplication of the digital data by the 
windowfunction, an interpolating processor which 
applies an interpolation processing by tertiary 
function approximation to the result of the fast 
Fourier transform processing and a display unit 13 
displaying the result of the interpolation processing. 
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